YITP-07-07, BARI-TH/07-563 



o 
o 

(N 



in 



BCS/BEC crossover in Quark Matter 
and Evolution of its Static and Dynamic properties 

- from the atomic unitary gas to color superconductivity - 

Hiroaki Abuki 1 ' 2 ^ 

1 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
2 I.N.F.N., Sezione di Bari, 1-70126 Ban, Italia 
(Dated: February 2, 2008) 

We study the evolution of dynamic properties of the BCS/BEC (Bose-Einstein Condensate) 
crossover in a relativistic superfluid as well as its thermodynamics. We put particular focus on 
the change in the soft mode dynamics throughout the crossover, and find that three different effec- 
tive theories describe it; these are, the time-dependent Ginzburg-Landau (TDGL) theory in the BCS 
regime, the Gross-Pitaevskii (CP) theory in the BEC regime, and the relativistic Gross-Pitaevskii 
(RGP) equation in the relativistic BEC (RBEC) regime. Based on these effective theories, we discuss 
how the physical nature of soft mode changes in the crossover. We also discuss some fluid-dynamic 
aspects of the crossover using these effective theories with particular focus on the shear viscosity. 
In addition to the study of soft modes, we show that the "quantum fluctuation" is present in the 
relativistic fermion system, which is in contrast to the usual Nozieres-Schmit-Rink (NSR) theory. 
We clarify the physical meaning of the quantum fluctuation, and find that it drastically increases 
the critical temperature in the weak coupling BCS regime. 

PACS numbers: 12.38.-t, 25.75.Nq 
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I. INTRODUCTION 



Quantum ChromoDynamics (QCD) is expected to ex- 
hibit a surprisingly rich phase structure at finite temper- 
ature and/or density. The lattice calculations strongly 
support the conjecture that the nuclear matter under- 
goes a phase transition to the color-dcconfincd quark- 
gluon plasma (QGP) phase at some critical temperature 
HLH. The QGP is a longstanding theoretical issue since 
the discovery of the asymptotic freedom of QCD 0] and 
is now being searched experimentally in the Relativistic 
Heavy Ion Collider (RHIC) [|. 

The QGP was originally thought as the weakly cou- 
pled plasma and its early conjectured signature [5(] was 
based on this observation, i.e., the strong suppression 
of heavy qq bound states above T c . However, the lattice 
calculations of mesonic spectral functions using the Max- 
imum Entropy Method (MEM) show that the mesonic 
bound states persist well above the critical temperature 
Q. Moreover, the RHIC data of the collective flow to- 
gether with the theoretical analyses using the parton cas- 
cade simulation and the hydrodynamics |U suggests 
the strongly correlated plasma where the partonic cross 
section is almost 50 times larger than its perturbative 
estimate. The recent lattice calculation of the shear vis- 
cosity to entropy ratio Q in the purely gluonic plasma 
also shows the value which is a factor 2-3 smaller than 
its perturbative estimate [To| ; the lattice value is rather 
close to the minimum bound (1/in) speculated in [llj 
using the string theory together with the AdS/CFT cor- 
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respondence [121 ]. 

These theoretical and experimental results show that 
the characteristics of the QGP are, (i) the existence of the 
quasi-bound states above T c and (ii) an almost perfect 
fluidity with small sound-attenuation length (l3l |. Some 
theoretical efforts have been devoted to understanding 
such a strongly coupled QGP 0, EI EE E3 in which 
there exist two temperatures, one for the deconfining 
transition or chiral-restoration (T c ) and the other for the 
vanishing of resonances (T* > T c ). The existence of the 
mesonic correlation above the chiral-restoration temper- 
ature T c was suggested earlier [l8| using the NJL model 
analysis [IE HE- 

The baryonic matter is also expected to be deconfincd 
to quark matter when it is compressed, which was con- 
jectured in plj or even earlier (22j. Moreover, a variety 
of color superconducting phases at high baryon density is 
suggested |23l . l24l [25| . It is now widely accepted that the 
ground state of quark matter is the Color-Flavor-Lockcd 
(CFL) phase [26| at extremely high density where the gap 
and critical temperature can be studied by the perturba- 
tive Dyson-Schwinger type equations [l?], [H, HE H3] ■ In 
contrast, which phase is realized at relatively low den- 
sity relevant to compact stars still remains a matter of 
debate because there is no general scheme. However, 
recent extensive studies based on effective models have 
revealed that there may arise a surprisin gly rich variety 
of non-BCS exotic phases [3l|, [H, [H, \34, l35l. l3rl l37l . [H| 
once the kincmatical constraints such as the neutrality 
with respect to the color and electric charges are incorpo- 
rated. Due to a relatively large strange quark mass, these 
kinematic effects plays an important role there. Such a 
stressed pairing is attracting broad interest not only from 
the QCD community but also from the audience of the 
atomic polarized Fermi gas [IE HE Sl|, |42j . 
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Another (and perhaps more direct) source of difficulty 
in investigating the pairing at low density is the presence 
of dynamical effects due to its strong coupling nature. As 
the system goes towards lower density, the gauge coupling 
grows and some strong coupling effects beyond the mean 
field approximation (MFA) come to play an important 
role as pointed out earlier by evaluating the Cooper pair 
correlation length within the MFA [43|,[44| . Some theoret- 
ical efforts to examine the strong coupling effects beyond 
the MFA were then made with a p articular focus on the 
precursory soft mode plU . liiil . l47l |. where non- vanishing 
diquark correlation and the pseudogap above the critical 
temperature are reported. Along with this line, the limit 
temperature T* where the resonances decouple from the 
spectrum is evaluated both for the chiral and diquark 
channels in p7| . 

Such existence of two temperature scales at strong cou- 
pling may be naturally understood in the scenario of the 
crossover from the BCS pairing to the Bose-Einstein con- 
densate (BEC), i.e., the BCS/BEC crossover [H, [H, H3] ■ 
The BCS/BEC crossover has been the longstanding the- 
oretical idea first discussed in the context of a theory 
of superconductivity in low concentration systems such 
as SrTiOs doped with Zr (48|. ft has also been intro- 
duced to understand the isotropic nature of elementally 
excitation in 3 He superfluid where the spin fluctuation 
results in strong attraction in 3 P channel leading to the 
non-s-wave pairing 14911 . Their work was extended to fi- 
nite temperature in [501 ] where the fluctuation about the 
MFA was inevitably included within the gaussian approx- 
imation. The basic concept of the BCS/BEC crossover is 
the following; in the weak coupling, the system exhibits 
the BCS superconductivity due to the attraction and the 
large density of state at the Fermi surface, while in the 
strong coupling, the composite bosons are formed at T* 
prior to their condensation to the bosonic zero mode at 
T c < T* . Although the symmetry breaking pattern is the 
same in both sides and therefore there is no sharp phase 
boundary in between, the mechanism of the condensation 
is completely different in the sense that the former has a 
dynamical origin while the latter has a rather kinematical 
origin; in the BEC side, the short range quantum effects 
are taken only into the structure of composite boson. 

So far, the BCS/BEC crossover has been widely dis- 
cussed in various contexts including not only the liquid 
3 He [49| or the high T c superconductivity [Hl|, but also 
the nuclear matter [53] , and the magnetically trapped al- 
kali atom systems j53[. Although the crossover scenario 
in color superconductivity was conjectured in the liter- 
ature [H, 0, H3 ! the first explicit investigation of this 
problem was given very recently in [55| where it is shown 
that the relativistic fermion system exhibits the two-step 
crossovers; one is the ordinary BCS/BEC crossover but 
with some unconventional behaviour of the critical tem- 
perature brought about by relativistic effects, and the 
other is the crossover from the BEC to the relativis- 
tic BEC (RBEC) [13 where the critical temperature in- 
creases up to the order of Fermi energy. Recently the ex- 



istence of the RBEC phase has been confirmed at T = 
in the Leggett's flamework and the evolution of collective 
modes are reported J5(| . Investigation of such a relativis- 
tic crossover is also performed in a boson-fermion model 
[58l . |59| which together with a chemical equilibrium con- 
dition enables us to describe the crossover thermodynam- 
ics with a boson mass or chemical potential controlled 
by hand. A possible importance of the diquark BEC 
in QCD phase diagram at low density is demonstrated 
in a model calculation [6(3|. The BCS/BEC crossover is 
also discussed in [6l[ with possible relevance to the chiral 
transition using the nonlinear sigma model. Moreover, in 
recent work jl6| , it is discussed with special emphasis on 
the chiral pseudogap phase above T c . 

As above, the BCS/BEC crossover is discussed in the 
QCD context to give an insight into the strong cou- 
pling nature of the quark-gluon plasma in either low or 
high density. However, how the transport and hydro- 
dynamic properties change throughout the crossover re- 
mains unclarified; this is interesting because it may pro- 
vide some unified view on the two apparently indepen- 
dent aspects of QGP, i.e., (i) the survival of the bound 
states above T c and (ii) the perfect liquidity. 

It is also worth mentioning that the hydrodynamic 
properties of the unitary Fermi gas have recently at- 
tracted considerable attention both experimentally (62l . 
IH, [H, [H, HI [H and theoretically [Hi, M, \M- The 
unitary Fermi gas is the intermediate of the BCS/BEC 
crossover, where the s-wave scattering length diverges. 
Such strongly coupled system can be created in the 
atomic traps with the external magnetic field fine tuned 
to the Feshbach resonance (7f[ . It provides us with a the- 
oretical challenge to describe such a unique many body 
system which is universal in the sense that there is no 
intrinsic dimensionful parameter except for temperature 
and density [72| . Understanding such system also may 
shed light on the physics of high T c superconductivity 
[Hlj . and possibly, the QGP. 

The aim of this paper is to understand how the 
transport properties of the relativistic fermion system, 
as well as its thermodynamics, change throughout the 
BCS/BEC crossover. In order to do this, we first formu- 
late the relativistic Nozieres-Schmitt-Rink (NSR) frame- 
work extending our earlier work (55| to fermions having 
SU(2)p flavor and SU(3) C color with a little refined regu- 
larization scheme. We then derive the effective theory for 
soft modes and study how this effective theory changes as 
the system moves from BCS to (R)BEC regime. Based 
on this effective theory, we will discuss some of fluid dy- 
namic aspects of the crossover with particular focus on 
the shear viscosity and its ratio to the entropy which 
serves as a measure of "perfect liquidity" . 

The rest of this paper is organized as follows: In Sec.HTl 
we formulate the Nozieres-Schmit-Rink theory for a rel- 
ativistic fermion system. It turns out that, in contrast to 
the nonrelativistic system, the relativistic fermion system 
has an additional source of fluctuation, i.e., the "quan- 
tum fluctuation" which we have just ignored in our previ- 
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ous study [5 51 ] . Our framework is not restricted to quark 
matter, and may be used for the other relativistic fermion 
systems such as a possible neutrino supcrfiuid trapped in 
compact stars at the early stage of their thermal evolu- 
tion [74j| . In Sec. IIII1 we discuss the static part of prob- 
lem, i.e., thermodynamics, the pair size, the spectrum, 
etc. In Sec. IIV1 we derive the effective theory for soft 
modes and study how it evolves with the crossover. In 
addition, we discuss the shear viscosity to entropy ratio 
based on the effective theory. In Sec. [V] we investigate 
the density dependence of the crossover. We also study 
how the "quantum fluctuation" affects the crossover. In 
Sec. I VII we make concluding remarks and perspectives. 



II. FORMULATION 



with G — N 4 " ] ^ 1 g 2 ■ Though we find the attraction in the 
scaler qq channel, we shall ignore this channel throughout 
this paper; consequently, our framework is not reliable in 
the vicinity of the chiral-restoration density. By intro- 
ducing a spinor doublet Q = (q,q) , we write 
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Here, we have defined P v — 70(^17)^70- Then the parti- 



tion function is 

Z = VQ exp 



' f 

cLt I cLxCe{Q) 



(5) 



In this section, we extend the Nozieres-Schmit-Rink 
theory so as to apply it to study relativistic fcrmions in- 
teracting via a point attraction. In Sec. Ill At we derive 
the thermodynamic potential up to the gaussian fluctu- 
ation about the MFA. In Sec. HTbI we introduce an in- 
tuitive representation of the thermodynamic potential in 
terms of the spectral density. In Sec. Ill C| we introduce 
the in-medium phase shift which enables us to handle 
easily the fluctuation effect to the number conservation. 
In Sec. lII D| we discuss the Thouless criterion and present 
how to renormalize it using the low energy scattering pa- 
rameter, the scattering length. In Sec. Ill El we give the 
renormalization scheme for the dynamic pair susceptibil- 
ity. In Sec. Ill Ft we end up with a couple of basic equa- 
tions, i.e., the number conservation, and the Thouless 
criterion. 



A. Nozieres— Schmit-Rink theory for a relativistic 
fermion system 

We start with the following four-Fermi model inspired 
by the one-gluon exchange (OGE) [55j ] . 

2 8 

C = q^ + ^-m)q- a -Y J ^T a 1 ^qqT a ^q. (1) 

a=l 

Here quark has N c = 3 colors and Nf = 2 flavors, and 
T a = A a /2 is the Gell-Mann matrix for color. At high 
density, we expect the pair formation in the J F = + , 
and the color- flavor triplet channel [23|, i.e.. 



(2) 



The indices {i, j, ■ ■ ■ } represent flavors and {a, 6, ■■•} in- 
dicate colors. Fierz transformation allows us to re-write 
the Lagrangian Eq. {1} as follows. 



C = q^ + ^-m)q+^j:l =1 [iq t P n q}[tqP n q t } 



, (3) 



where Ce = — £|t_>_, T is the Lagrangian density in the 
Euclid space, which is periodic in the imaginary time in- 
terval [0,/?]. Introducing the Hubbcrd-Stratnobich fields 



(t, a;) — 2 ( Q [ p nl^ 





Pr, 

P, 




Q 



(6) 



and integrating out quark fields, we obtain 

IAHt,*)! 2 



z = 



VAVA* exp 



x exp 



irdx 

-1 



G 



ilogDetf P ^ 5 -y 



S 



Fx.y 



Here the bare propagator is defined by 



S fI, v = + ^-m)S x . y , 



0-1 

°Fx,y 



75^^75- 



(7) 



(8) 



6 x . y is the delta function antipcriodic in the imaginary 
time r; in momentum space, 



Jx,y 



dp 



-iu) n (r :c -T y )+ip-(x-y) 



(9) 



with io n being the fermionic Matsubara-frequency. In- 
troducing the following notations for the Nambu-Gor'kov 
propagator and self energy, 



Fx,y 



-^x,y — 



s 



Fx,y 







0-1 

^Fx,y 







(10) 



we can write the partition function in the following form 
where the fluctuation contribution is factorized: 

ZQi, T) = e -^ T ) = Z (fi, T)Z fluc (/i, T), (11) 
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where Zfi uc is the part coming from the fluctuation: 

|A r; Hr,:E)| 2 



VAVA* exp 



drdx 



G 



x exp 



-log Det ( 5 x ,y - ^ SFx,z^z,y 



J2 Z is a shorthand notation of dr z J dz. The free par- 
tition function Z is e~^ n ° with f2 being the thermody- imation), we have 

I 



namic potential for free (massive) quarks: 



no 
v 



a=± 



-e pa /T 



(13) 

where e p ± = m 2 + p 2 =p p, and the vacuum fluctuation 
(12) which has no dependence on (/z,T) is suppressed. 

Up to quadratic order in (A, A*) (the gaussian approx- 



-/3n fluc 



II / dA v (in N ,P)dA*(in N ,P)exv fL- XlltT (iQ Nt p) \ \A v (m N ,P)\ 2 

AT D L \ / 



(14) 



where fi^v and P denote the bosonic Matsubara-frequency and momentum. The pair correlation function X/j.t(^tv, P) 
at one-loop level is defined by [HI, [T^ 1 



X*i,T(*fitf,P) = 2Tj] J-0^tT^S F {iw n +iQ. N ,q + P)S F {-iUn,-q) 



= -2 
-4 
+2 



(2tt)' 



l | m 2 + q-(q + P)y- f F (E q+P - fi) - f F (E q - f i) 

iQpj + 2/i — E q+ p — E q 



E q E q+P 



f i _ m 2 + q-(q + P) \ -f F {E q+P - p) + f F {E q + M ) 



J(2tt) 3 \ ± E q E q+P J in N + 2ii-E q+P + E q 

f dg / m 2 + q-(q + P) \ 1 - f F (E q+P + (j.) - f F (E q + 
7 (2tt) 3 V E q E q+P ) itt N + 2 l x + E q+P + E q 



(15) 



where S F (uj n , p) is the Fourier transform of the fermion 
propagator S Fx , y , i.e., the inverse of Eq. ©. Integrat- 
ing out the gaussian fluctuation leads to the following 
expression for thermodynamic potential. 



n0i,T) = n 0i ) T) + n flno (/i,T) ) 



(16) 



The dynamic pair susceptibility T(lu, P) can be ob- 
tained by the analytic continuation of pair correlation, 
Eq. (fT5|) , to the real w-axis: 



t;I f (uj,p) = ^-x,,t(u + iS 1 p). 



G 



(18) 



with the gaussian fluctuation being defined by 



fiuc 



dP 
dP 



log 



1 

G 



Xfi,T(itt N ,P) 



(2^) ; 



log 



--Xo(ifi,P) 



where ds — 
son and xo 



N C (N C 



(17) 

is the number of "flavors" of bo- 
Xfj,=a,T=a] we have subtracted the "vac- 
uum" contribution at fi = T = 0, which in general leads 
a G-dependent fermion mass and wavefunction renor- 
malization through a Hartree and higher order terms 
[75l | . Some notes are in order here, (i) The fluctu- 
ation is symmetric under fi — (i as it should be. 
This is guaranteed by the charge conjugation property 
X-p.t(-«^w,P)| = Xm,t(^JV,P)- (ii) The above for- 
mula is analytic above some T c below which the system 
undergoes the transition to pairing phase. 



B. Thermodynamic potential in terms of spectral 
density 

We introduce a spectral density at some coupling 
strength G = Q by 



i/g -xM^n,p) 



ir ifl^ — u 



(19) 



It is clear that the spectral density is related to the imag- 
inary part of the dynamic pair susceptibility at Q. 



pj* (o,,P)=ImrS (w,P) 



(20) 



To express the thermodynamic potential in terms of the 
spectral density, we differentiate and integrate Eq. (fT7| 
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with respect to G. We obtain 
r G 

^fiuc = d B 



dg 
Q 2 



, 7T (27r) 3 

-(T = p = part) 



du dP T ^ upl >T {u : P) 



N 



n 2 



(21) 

Because of the identity p_^ T (w, P) = — p 5 T (— u, P), 
the thermodynamic potential still possesses the charge 
conjugation symmetry. By performing the Matsubara- 
summation, we obtain 



o G 

G 



du dP ~ 



dg 

2 '-< 

dg 



du> dP e(u>) . a , 

>Ap s (u,,P). 



Q 2 J_ x 7T (2tt) 3 2 

(22) 

where f B (w) = f B {u) + 6{-u) with f B {u>) = l/(ef u - 1) 
being the bose distribution function. We have defined 
the in-medium spectral shift Ap^ T (to, P) = p^ T (uj,P) — 

Pq(u>,P) with Pq(lo,P) = Pq (lo,P) being the spectral 
density at p. = T = 0. We denote the first term by 
^NSRi and the second term by Q q a hereafter. The total 
fluctuation is then the sum of these two pieces. 



n 



NSR 



qfl • 



(23) 



where 



^NSR 



and 



G 



dg 



(24) 



f2 q fl = — dl 



G dg 
g 2 



duj dP e(oj) 
,7 (27r) 3_ 2~ 



^ r (w,P). 



(25) 

f^NSR is the Nozieres-Schmit-Rink correction to the ther- 
modynamic potential and roughly corresponds to the 
thermal fluctuation, while f2 q fi can be regarded as the 
quantum fluctuation which is ignored in the nonrelativis- 
tic Nozicrcs-Schmit-Rink theory. If we take the limit 
T — > 0, the former goes to zero, but the latter remains: 

-Mo G f/:!^ [p%{u,p)-p1^p)\ . 

(26) 

fiqfl may play an important role at low temperatures, 
but we shall ignore this contribution for a while because 
(i) the fluctuation itself is less significant in the weak 
coupling (low T c ) regime, and (ii) the charge conjugation 
symmetry is maintained even if il q fi is ignored. We will 
come back to this problem in Sec. [V] where we examine 
how large this "quantum correction" is. 



Thermodynamic potential in terms of phase 
shift 



It is sometimes more convenient to express the ther- 
modynamic potential in terms of the in-medium phase 



shift rather than the spectral density. It may be defined 
by the integration of the spectral density: 



G dg 
) Q 



%(% iT (LJ t P) = -log 



%-X»,T(w + iS,P) 
5 - Xti.riu - iS,P) 



<Wr(w,P), 



i.e., the argument of the dynamic pair susceptibility 



(27) 



X^.t(uj,P)\ 



(28) 



The phase shift possesses the charge conjugation 
S-ii,T(—v,P) = — S(i,T(u, P )■ Using this phase shift, 
we can express JInsr an d £^ q fl as 



NSR 



n 



qfl 



-oo * (2tt) 3 
dto dP e(oS) 
,V(27r) 3 2 



(29) 



A8^ T {u),P), 



where the in-medium shift of the phase shift is defined 
by A5^ t {oj, P) = 5^ T (u, P) - 8 (u, P) with S (u, P) = 
6qo(u>,P) being the phase shift in vacuum, ^nsr is ex ~ 
actly of the same form as that derived by the Nozieres- 
Schmitt-Rink (H. 



D. Thouless criterion and the renormalization of 
the coupling 

In the weak coupling regime, the Thouless condition 
determines the critical temperature. This condition en- 
sures the divergence of the long wavelength limit of the 
dynamic pair susceptibility at the critical temperature 
for a given chemical potential p: 

r ;,T c (0,O) = i- XM , Tc (0,O) = 0. (30) 



Xju,t(0, 0) can be calculated as 

Wo, oh 2 



dq tanh ^fp^- 



3 E q 



P 



+ {p^-p). (31) 



which is quadratically divergent. Accordingly, the critical 
temperature should depend on the cutoff as 



T c = A/(m/A,/i/A,GA 2 ), 



(32) 



where / is some function. We can reduce the cutoff de- 
pendence via the partial renormalization of the coupling 
using the low energy information about the scattering 
T-matrix. We here consider the fermion-fermion scatter- 
ing (/(l) + f(2) -> /(3) + /(4)) with (1, 2, 3, 4) labeling 
the momentum, color, flavor and hclicity of each quark, 
as ("1" = (p,a,i,hi),"2" = (-P,b,j,h 2 )), ("3" = 



6 



(k, c, k, /13 



-k,d,l,hi)). Then the on-shell T- 



matrix can be parameterized as 

T(12 - 34) = T(p, fe) (r 12 r 34 - (3 «-> 4)) , (33) 
where T 12 is defined by 



£gb gjj (T h 1 h 2 



(34) 



<r 3 is the Pauli matrix. T(p, k) can be evaluated with the 
Lippman-Schwingcr equation as 

-G 



T(p,k) 



= -r (2E p ,0). (35) 



Here, To is the dynamic pair susceptibility, Eq. (|18p . at 
/i = T = 0, i.e., T (cj, P) = r ,o(w, P). Using the defini- 
tion of the phase shift function introduced in the previous 
section, we find the following expression for the scatter- 
ing amplitude f(p,k) = — ^T(p,k) with (m = to/2) 
being the reduced mass 

f(p,k) = -^-\\T (2E p ,0W 5 ^ o l (36) 

At sufficiently low energy p = k <C to, one can make 
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use of 2E p = 2m + ^ in the above formula, and the 
result should be matched with the general form of the 
low-energy expansion of scattering amplitude, i.e., 



/(P,fc) = - 



pcotS— ip — + ir e p 2 — ip 1 (^^) 



where a s is the s-wave scattering length and r e is the 
effective range. We find the scattering length is quadrat- 
ically divergent as 



to 1 
= G 



Ana 



1 



dq 

(2tt) 3 \E q - to ' E q 



(38) 

From now, we use the renormalizcd coupling Gr instead 
of a s itself for notational simplicity; that is defined by 



1 



Eq. 



G R 4na s 
can be written as 



1 

Gr 



1 

G 



Xo (2to,0). 



(39) 



(40) 



Using this renormalizcd coupling, the gap equation can 
be casted into the following form: 



dq ( tanh 



1 



(2tt) 3 V E q - /x E q -m ( 41 ) 
+ (yU, to — > — /!, —to) 



1 



It can be easily seen that this integral still has a weak 
logarithmic divergence. This is in contrast to the non- 
relativistic case where the UV divergence can be com- 
pletely taken away because of the quadratic momentum 



dependence of the single fcrmion excitation, p 2 /2m, in 
the energy denominator. T c is now parameterized by 



T c = mf R (ij/m,GRm 2 ;\n(A/m)). 



(42) 



Jr is again some unknown function. The apparent cutoff 
dependence seems to be smaller than the Eq. (pi2l . 
Eq. (|4Tj) can also be formally written as 



1 

Gr 



xj&(o,o) 



0, 



(43) 



with the renormalized pair correlation function 



P) = xM^ P) - Xo(2to, 0). (44) 

We next discuss the bound state equation in vacuum 
fi = T = in terms of the renormalizcd coupling. The 
bound state or the resonance pole in two fermion scat- 
tering can be determined by 



J-- X £e> )0) 



0, 



(45) 



where x? en (w, P) = Xo^O, P)- Because x? en (2m, 0) = 
0, if —1/Gr > then the resonance pole is located at 
I a; I > 2m. Otherwise it is located at \u>\ < 2m, which 
corresponds to the stable bound state because Imx(M < 
2m, 0) = 0. Then the critical coupling Gr = Go for the 
zero binding is given by the condition 



1 



r< Xo 



Ren (2TO,0) = -— = 0, 



f 



(46) 



i.e., l/a s — which is sometimes called as the uni- 
tary limit. If the coupling is above this critical coupling 
(1/Gr > 0) and if the bound state is sufficiently loosely 
bound so that the pole u> satisfies < 4to 2 — uJ 2 <C to 2 , 
then the bound state pole can be approximated model- 
independenty by the unitary condition; that is 



2to 



1 



2to 



16^ 2 
m 3 G 2 R ' 



(47) 



The binding energy is given as usual by 1/ma 2 . In our 
relativistic case, we have another typical coupling l/G c 
which is stronger than the unitary coupling: When the 
coupling approaches this critical point 1/Gr — > 1/G C — 0, 
the bound state becomes massless. This critical coupling 
satisfies the following condition, 



G c 



x? en (o,o) 



lx? en (o,o)|. 



(48) 



When the attraction is increased beyond this G c , the 
vacuum becomes unstable against the formation of qq- 
condensate which leads to the Majorana mass gap in the 
single quark excitation even in vacuum. 

It is also noteworthy that because of the identity 



X Re > + 2 M ,0), 



(49) 
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the Thoulcss condition (-1/G R - X^t(°>°) = °) 
rather determines the chemical potential \i in the region 
-1/Gr < and 2(m-/i) > T where the Pauli-blocking 
gives only a minor effect (x^(0, 0) - xffi(0, 0)). Such 
conditions are actually realized in the nonrelativistic 
BEC regime, where fi should be about one half of bound 
state energy at rest, but with small T-dependent entropic 
correction [7a ]. 

Before closing this section, let us briefly discuss the 
to — > limit. Although we have introduced the renor- 
malized coupling Gr via a s , the low energy informa- 
tion about f(p, k) assuming p <C to ^ 0, if we extend 
the definition of Gr by Eq. (|4Tj) . the succeeding discus- 
sion can be applied even to the m = case. By defin- 
ing an appropriate function f' R , Eq. ([42]) can be casted 
into T c = /i/^(m//i, Gr{j, 2 , ln(A//i)) in which it is rather 
clear that it has a smooth massless limit. Also in this 
limit, 1/Go = 1/G C = as is clear from their definitions, 
Eq. gSl) and Eq. ((48]) combined with x^ on (0,0) = 0. 

E. Renormalization of the dynamic pair 
susceptibility 

In this section, we polish up the argument in the pre- 
vious section. By doing this, it turns out that we only 



need a regularization of the vacuum (and also real) part 
of the dynamic pair susceptibility. 

First we notice that the imaginary part of r~]p(w, P) is 
finite while the real part is divergent. Thus, we only look 
at the real part Re r~y(o;, P) leaving the explicit formula 
of the imaginary part in the Appendix. The real part of 
the pair dynamic susceptibility can be decomposed into 
two part: 



ReT-^u, P) = Rer- V + 2/x, P) - Rex^(w, P), 

(50) 

where we have defined the "pair susceptibility" in the 
vacuum (/t = T = 0) as 



r -1 (^) = y;-Xo(z,P), (51) 



with xo being Xh=o,t=o- The explicit forms of \o an d 
X™t can be- read from Eq. (fT"5|) . 



Re X o{z,P) = -2 



= +4 



dq 

dq 

(2^)3 

dq 
(2nf 



1 - 



m z + q-(q + P) 
E q E q+P 

m 2 + q ■ (q + P) 
E q E q+P 

to 2 + q ■ (q + P) 

EqE q+ p 



V 



V 



E q +p — E q Z + E q+ p 



E n 



Vf F (E q -fx) 



Vf F (E q + p) 



u + 2/i 

Vf F {E q +fx) 



E q+P - E q oj + 2/i + E q+P + E q 



Vf F (E q -p) 



2/i — E q+P + E q 



2/i — E q + E q+P 



(52) 



From these formulas, we can see that x™t ( w i -P) i s finite 
while the vacuum part xo( z yP) is quadratically diver- 
gent. In the same way as in the previous section, we can 
reduce the divergence in T^ 1 by the renormalization of 
the attractive coupling (see Eq. I|38p ). 



Rer- 1 (o. + 2/ / ,P) = I- Xo ( w + 2/i,P) 

- -^--X? e > + 2M,P) ^ 

(jR 

= Rer-jJ eii ( W + 2M,P), 



with xo en (w + 2/i, P) = xo(u + 2/i, P) - Xo(2m, 0). The 
momentum integral in the right hand side is again loga- 
rithmically divergent. We arrived at the following rcnor- 



malized pair susceptibility: 

[r*-(a,,P)r a = Rcr- R 1 cn (^ + 2/i,P) 
-Rex™f (w,P) 
-iImx MiT (w + i^ P). 



(54) 



We note here again that only the first term has a loga- 
rithmic divergence. The renormalizcd Thouless condition 
which we have derived in the previous section can be ex- 
pressed in the compact form: 



[r^^o)]- 1 



o. 



(55) 



This is exactly the same as Eq. (TT3]) if the matter part of 
the integral in Eq. (|13[) . the piece being proportional to 
fp, is evaluated without a cutoff A. 
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F. Number conservation 

Using the renormalized dynamic pair susceptibility dis- 
cussed in the previous section, the in-medium spectral 
density and phase shift can be defined by 



S%p(w,P) = Arg[r^( W ,P)]. 



(56) 



The total density AT tot (//, T) is given by the derivative of 
the thermodynamic potential with respect to /i: 



N tot (^T) = iV MF (/i,T) 
+d B 



du dP - . dSffifaP) 
— f B {uj) — 



7T (2?r) 3 



(57) 

where Nmf{h, T) = — ®Q°lti±Ll [ s the quark number den- 
sity from free quarks. This contribution can be further 
decomposed into two parts, i.e., the contribution from 
quarks, and that from antiquarks, i.e., 

Nmf(^T) = N q (ji,T) - N q (n,T), (58) 

with 



lj0ff F (E q -fx), 



N q {n,T) = 2N c N f 



dq 



(59) 



The second term in Eq. Q57[) represents the quark number 
density from pair correlation, which we denote by 



'dej dP - dS^u^P) 



7T (2tt) 3 



(60) 

In order to give an intuitive interpretation of this for- 
mula, we extract some specific contributions from this 
expression. We first note the identity 



d^ c T >,P) 



dRcfr^^P)- 1 ] 



dim rr^fcj.p)- 1 ! 

-Re Fffi (w,P) 1 ^' TV ; J 



(61) 

We focus on the first term because the second term plays 
only minor role as long as (fi, T) satisfies the Thouless 
condition. The number density from the first part is 



(62) 

If the attractive coupling is strong enough (1/Gr ^ 0), 
we have a stable bound boson (antiboson) pole in the 
spectral function. 



po lc p ai , = Z£ T *(« + 2M-2#£) 



-Zf8^ + 2 l i + Ef p ), 



(63) 



where Z^ T , E^p (Z^ T ', E^ p ) arc the wavefunction 
renormalization and the bound state energy for boson 
(antiboson). E^p and E^p are defined by the following 
bound state pole conditions: 



r^T n (^'p-2 M ,p)- 1 = o, 

I^(-^-2/i,P)- 1 = 0. 



(64) 



Imaginary parts should vanish for stable bound states. 
The wavefunction rcnormalizations are calculated by 



Tt*,T 



-2j.i 



(65) 



-2/i 



Substituting these expressions to Eq. (f6"2"| . we find the 
following bound state contributions to Nnsr(^,T) 



r dP 


urj BP 


1 (2^) 3 


dfi 



r dp 


2H 


BP 


I (2tt)3 




dfi 



(66) 



f B {Ef p +2fi). 



We denote the first term by N B {fi,T) and the second 
term by —N B (fj,,T). Note that the "effective" quark 
number charge of (anti)boson slightly deviates from (— )2, 
which is caused by the Pauli-blocking effect. The remain- 
ing contribution to A^nsr(m, T) comes from unstable pair 
correlation such as the ^-continuum excitation or the 
Landau damping. We define these contribution by 



N un (fi,T) 



N mR (n,T) - N B (ji,T) + N B (n,T). (67) 



The chemical potential fi c may be determined by the 
quark number conservation: 



k 3 

Ivp 



N c N f ^ = N tot (fi c ,T) 



AW(At c , T) + N NSR (fi c , T), 
(68) 

where the Fermi momentum kp is introduced as a param- 
eter which controls quark number density; we fix the total 

k 3 

charge of the system as Ntot = N c Nf^s- In the weak 



3n 2 

coupling, this equation actually determines the chemical 
potential fi c ~ Ep = \/k'p + m 2 , while in the strong cou- 
pling region, it rather determines T c because the Thouless 
condition determines fi c . 

In the following, we derive some analytic approxima- 
tions for T c in the strong coupling. As noted, the Thou- 
less condition in the strong coupling is essentially the 
bound state (Bethe-Salpeter) equation in vacuum where 
2fi c plays a role of boson mass M B = 2fi c < 2m. Then 
the fermionic contribution to the number density gives 
only a minor contribution and instead the bound state 
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Weak rounline RCS (C1~ L <S 01 


udou v y^L ± tiling 


RRC CO <5t (7~ x <st (7~ x 1 


(A) Thouless criterion 


T c ~E F e 2fc Fi»= 


(Mc,Tc) 


2u c ~ 2m 


(B) Number equation 


fj, c ~ _Bf 




J> - . k F ( 2 V /3 
c 7rV 3 m ^3(iV -l)C2/3 / 


Physical gap at T = 


2A ~ 27re^ 7E (= 3.53)T C 


~2\/A ( 2 ) + (m-/i c ) 2 


2m — 2/i c 



TABLE I: The roles of two basic equations. Ef = \Jk F + m 2 is the Fermi energy, and the s-wave scattering length is related 
to the renormalized coupling by a a — tuGrI^k. 



contributions dominate the total density. Therefore we 
approximate the number equation Eq. (|68[) as 



2d r 



dp 



(2tt)= 



f B (E%f - 2/i c ) - f B (E%f + 2n, c ) 



(69) 

We further approximate this for the following two cases, 
i.e., (1) nonrelativistic case and (2) relativistic case. 

(1) Nonrelativistic BEC; this means that the boson 
mass 2/x c is much larger than the boson's kinetic energy 
or T c , i.e. 

2 Mc » N^/Aftc. (70) 
This may allow us to approximate boson's dispersion as 



2jJL c 



p2 

4^ c 



(71) 



Then the boson pole in Eq. (|6"3j) dominates the integral 
because of the bose-enhancement factor. We have 



2d f 



7T 



3/2 



C(3/2). 



(72) 



Therefore, if the temperature becomes smaller than the 
critical value 



t ^ t nr — 

1 ^ J BEC — 



N f 



1 



tt 1 ^ ViV c -13C(3/2) 



2/3 



(73) 



the macroscopic fraction has to be condensed in the P = 
mode to maintain the number conservation; that is 



the classical BEC for the nonrelativistic bosc gas. The 
critical temperature is of order N^/Afi c (<c 2u. c ) which 
also justifies our approximation Eq. (|7ip . 

(2) Relativistic BEC; there is another situation which 
was first considered in [53]; that is the case 



2fi, 



« < 7 t 3 . 



(74) 



In this case, we have to take care both of the boson and 
antiboson contributions in Eq. (|69[) . We may approxi- 
mate the boson and antiboson dispersion by \J Au. 2 . + P 2 . 
Then, we obtain the following approximate formula in 
the totally same way as (53|, 



2/i c T 2 



2d B 3 • <™ 

Then the critical temperature to the relativistic BEC 
(RBEC) is given by 



1 BEC — 



Nf 



TT V Mc 2(JV C -1)' 



(76) 



The critical temperature is of order (TVtot/2/^c) 1 ^ 2 which 
is much larger than the boson mass 2\x c . 

The crossover regime: If the coupling is in the in- 
termediate range such that the system is in the crossover 
regime between the BCS and the BEC, the coupled set of 
the Thouless condition and the number equation deter- 
mines (fi c ,T c ). Thus, our basic equations to determine 
(/i c , T c ) in the whole region of the attractive coupling are 



(A) Thouless condition: = Re [r^ n Tc (0, 0)" 1 ] , 

(B) Number conservation: N c N f ^ = AT MF (/i c , T c ) + J -^J^IbH 1 ^ ' ■ 



We summarized the roles of the Thouless criterion and 
the number equation in the BCS and BEC regimes in 
TABLE. [T] It may be worth noting that the physical 
gap at T = for the charge neutral excitation is always 
handled by the Thouless condition, i.e., by the low energy 



(long wavelength) limit of the dynamic pair susceptibility. 
We will see that, in the crossover regime, or in the ultra 
strong coupling regime, these two equations are strongly 
coupled, giving rise to a complex behaviour of (/i c ,T c ) as 
a function of the renormalized coupling Gr. 



III. THERMODYNAMICS OF THE 
BCS-TO-BEC CROSSOVER 

Here, we discuss how the thermodynamic character of 
the system changes throughout the crossover. In the nu- 
merical calculation, we set iV c = 3, to/A = 0.2 and fix 
the fermion number density at kp = 0.2to. 

A. The BCS, BEC, and RBEC phases 

FIG. HJa) shows how the critical temperature and 
chemical potential (T c and fi c ) change along with the 
attractive coupling G c /Gr; we use G c , the critical cou- 
pling defined by Eq. (|48|) . as a normalization. Figure 
(b) shows the quark number contents. We can clearly 
see that there are three distinct regions, which we have 
called the BCS, BEC and RBEC phases in our previous 
paper [55| . The critical temperature in the BCS region 
is well-approximated by the mean field result (T^; thin 
grey line) where T ( ^ IF and /i^ F are determined without 
the fluctuation effect to the number conservation, i.e., 
the result in which the last term of (B) of Eq. ([77)) is 
ignored. In this regime, the weak coupling universal re- 
lation between the zero temperature gap Ao and T c , i.e., 
T c ~ 0.567Ao is well realized as we can see in FIG. [2j 

As G c /Gr grows, T c gradually deviates from the mean 
field result and for G c /Gr 0.07 the bound state for- 
mation takes place and the system goes into the BEC 
regime. The bound state formation in medium (kp =/= 0) 
always takes place at a stronger coupling than G c /Gr = 
0, the unitary coupling. This is because the Bethe- 
Salpeter kernel is smeared the Pauli-blocking effect at 
finite density. FIG. [3] also shows this fact. The bold 
line shows the half of boson mass Mb(^ c ,T c ) /2 = /i c 
while the thin line represents the half of antiboson mass 
Mg(fi c , T c )/2 as a function of Gr. We can clearly see 
that the antiboson forms prior to the formation of boson. 
The antiboson forms almost G c /Gr ~ where the bound 
state forms in vacuum, while boson does not form up to a 
little stronger attraction G c /Gr = 0.07; this means that 
the Pauli-blocking by on-shell fermions actually prevents 
the formation of boson. 

In the BEC regime, the growth of T c as a function of 
G c /Gr is suppressed because the increase of attraction is 
mainly used to reduce the in-medium bound state mass. 
Because the in-mcdium boson mass changes, temperature 
does not saturate to the ideal BEC temperature of boson 
with mass Mr = 2m which is indicated by the arrow in 
the figure; this is in contrast to the nonrelativistic calcu- 
lations [IjMIizI- 

T c in this region is well approximated 
by the nonrelativistic ideal BEC temperature (T^ c ) of 
the boson with mass 2/i c , 

kp 1 / 1 \ 2 ^ 3 
t bec = k F ttt ^7 — • (78) 

fl c 7T ' V 3 C2/3/ 

This is shown by the thin orange line in figure (a). It 
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-0.5 0.0 0.5 1.0 1.5 



FIG. 1: (a) The critical temperature (chemical potential) 
T c (fi c ) as a function of attractive coupling G c /Gr. The dis- 
sociation temperature of pre-formed pair is also depicted by 
dashed line, (b) The quark number fractions. 



is remarkable that this expression has no explicit depen- 
dence on Gr; there is only the implicit dependence on 
Gr thorough 2^ c , i.e., the in-medium boson mass. This 
fact expresses that the condensation in this regime has 
a rather kincmatical origin. The factor N c — 1 in the 
above expression also indicates the condensation in this 
regime depends on the kinematical degrees of freedom; 
the factor comes from the fact there are N c fermions in 
the system while we have ds = N C (N C — l)/2 bosons 
belonging to the 3 C representation of SU(3). This is in 
contrast to F in the BCS region where we find no para- 
metric dependence on N c . It is also noteworthy that the 
universal relation T c / Aq = 0.567 is significantly violated 
as seen in FIG. [2 T c /Ao grows up to ~ 10 5 . However, 
the physical gap in the single quark excitation at T = 0, 
y/ Aq + 6(m — /i c )(TO — /i c ) 2 , is always comparable to T c . 

As the attraction, G c /Gr, is further increased beyond 
the unity, the system goes into the RBEC phase where 
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FIG. 2: The zero temperature gap parameter Aq, and T c of G^, 1 . Sf = iV c JV/( Mc 3 ° + 



as a function of G fl . The mean field gap parameter Ao is 
calculated at (/x = /u c , T = 0). The physical gap in the single 
quark excitation at T = 0, \J Aq + max.(m — fi c , 0) 2 , is also 
shown by the bold red line. 



FIG. 4: The entropy par quark number charge as a function 

is the entropy for free 
1)— tt*- is the Stephan- 



45 

massless quark gas and Sb = N C (N C 
Boltzmann entropy for N C (N C — l)/2- "flavored" boson gas 
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FIG. 3: The in-medium boson and antiboson masses 
(Ms,Mg) as a function of attraction . The cross indi- 
cates the point G c /Gr = 0.07 where the bound boson forms, 
i.e., Afs(/i c ,T c ) = 2/j.c = 2m holds at this point. 



T c exceeds the in-medium boson mass 2^ c . The critical 
temperature in this regime is well approximated by the 
ideal BEC temperature T^ c of a relativistic bose gas: 



rpKL _ ™F 

J BEC " 




(79) 



critical temperature, there are a lot of bosons and anti- 
bosons as shown in Fig.[T](b). We notice that unlike the 
ideal BEC of the elementary boson gas (57[, there are a 
lot of (anti)fcrmions as well as (anti)bosons in the RBEC 
regime. This is the characteristic feature of the com- 
posite boson system, where there exists the competition 
between free energy and entropy; two fermion state is 
more favorable than one boson state in terms of entropy, 
but is less favorable in terms of free energy. This con- 
sists one reason why we find a relatively large deviation 
between the real T c and T^g C in the RBEC regime. 

To confirm the above-discussed physical picture, we 
show in FIG.[U how the total entropy behaves going from 
the BCS to the RBEC. We can see that the entropy in the 
BCS regime is suppressed from the free massless fermion 
result because of the pairing correlation and finite mass. 
When the system goes into the BEC region, the entropy 
takes a nearly constant value. This is because the in- 
crease of the coupling mainly affects the internal struc- 
ture (the binding) of the boson and the kinetic degrees 
of freedom controlling T c does not change. When the 
system goes into the RBEC phase, the entropy rapidly 
increases because there appear a lot of kinetic degrees of 
freedom, {q, q, B, B}. Since plenty of {q, q) are present in 
addition to (B,B), the Stephan-Boltzmann entropy for 
the ideal bose gas underestimates the total entropy. 



B. Character change of Cooper pair and 
Deformation of Fermi surface 



To obtain a more intuitive insight into the three re- 
Because the in-medium boson mass is smaller than the gions, we here study the character change of Cooper pair 
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wavefunction, the spectral density, and fcrmion occupa- 
tion numbers. In FIG. [5J we show the spectral density 
[(a); left panel], the Cooper pair wavefunction [(b); mid- 
dle panel], and the (anti)fermion occupancy [(c); right 
panel]. From top to bottom, the renormalized coupling 
grows as G C /G R = -0.35 (BCS), G C /G R = 0.5 (BEC), 
and G C /G R = 1.35 (RBEC). 

We first discuss how the spectral density evolves 
from the weak coupling BCS to the strong coupling 
RBEC phase. Figure (al) shows the spectral density 
rri 2 p^ eu T (w, 0) in the BCS phase. We can see that the 
Thouless "zero" of the inverse susceptibility at to = 
is located above the 2-particle "decay" (and 2-hole "ab- 
sorption") continuum threshold, i.e., > 2m — 2\x c . This 
means that the pair fluctuation with finite momentum 
in the BCS state decays into two quarks. We confirm 
this later by constructing the low energy effective the- 
ory for the fluctuating pair field. FIG. [5] (bl) shows 
the spectral density in the BEC state. In contrast to 
the BCS case, the Thouless singularity at lu = is re- 
alized as a bound state (boson) isolated pole. This is 
because u> = is located below the 2-particle continuum, 
i.e., < 2m — 2fi c . There exists another singularity in 
the region to + 2/i < corresponding to the in-medium 
antiboson pole. We note, however, that the thermody- 
namic quantity should be expressed by the integral of 
the combination x p^^(uj,P), so the antiboson 

pole gives only an exponentially suppressed contribution 
to the thermodynamic quantities like entropy, heat ca- 
pacity, etc. If we go higher temperature T > T c with 
fixing jjb = p c , then the boson pole gradually shifts to 
higher u) because the bound state dissociates thermally 
(see the thin orange line in the figure). As we increase 
the temperature, the pole eventually gets absorbed into 
the 2-particle continuum at the dissociation temperature 



T 



The spectral density at T 



is de- 



picted in the figure by the thin red line. Interestingly 
enough, the antiboson pole still survives at this tempera- 
ture; this is attributed to the fact that the Pauli-blocking 
effect is less significant in the antiboson sector than in 
the boson sector. Finally, we show the spectral density 
m 2 /9^ n Tc (w, P = m) in the RBEC phase in FIG. M (cl). 
Because fi c <C T c in this regime, the spectral function is 
almost antisymmetric with respect to lo «-> —u>. In the 
region — P < u + 2p < P, the decay of the pair fluc- 
tuating mode by absorbing thermally excited antiquarks 
(the Landau damping process) is kincmatically allowed. 
Because of this, the spectral density takes non-zero value 
in this region. In the RBEC phase, the temperature is so 
high (T c ^> He) that this Landau damping process gives 
a significant contribution to the unstable content N un in 
^nsr(Mci^c)- This is one of the characteristic features 
in the relativistic fermion system. 

Next we study how the internal structure of the Cooper 
pair changes throughout the crossover. Let us first begin 
with the definition of the Cooper pair "wavefunction" . 
We may define it from the expression of the correlation 



functions for T <T C [73]: 

(qf(0,x)q b j (-i0-,y)) 
= +(r 2 X 2 )^ J j0^e^-y^ + (q)A 75 CA + (q) {gQ) 
+(r 2 A 2 )?// I ^e^-(-«)^_(g)A 75 CA_(g), 

where the quark-quark pair wavefunction (f> + and the 
antiquark-antiquark pair wavefunction arc defined by 



4>±(q) = tanh 



'<! + 



1 



2T J 2e q± ' 



(81) 



e q ± = y/ (E q =F p) 2 + A 2 is the quasi-(anti)quark en- 
ergy. These wavefunction in q-space can be Fourier- 
transformed to the real space: 



<j)±(r) 



1 



27r 2 r 



(82) 



Unlike the nonrelativistic situation, this integral does not 
converge even for r / 0; the integrand oscillates with 
sin(gr) as q — > oo. Evaluating the oscillating integral 
by the sharp cutoff A is inappropriate, and therefore we 
take the following prescription to regulate this integral. 
By integrating by parts, we obtain 



i±(r) 



l 

" 2-rr 2 r 2 



q cos(qr) 



tanh 



dqcoB(cpr)& 



2^tanh(^; 



(83) 

The second term is conversing while there remains a os- 
cillating uncertainty in the first term. However, if we 
adopted the smooth cutoff scheme, the first term should 
vanish. Thus we may define the second term of Eq. ([83]) 
as the regularized Cooper pair wavefunction: 



0±(r) 



i 

~2tt 2 t 2 



r o °°dgcos(gr)- 



tt— tanh 

2e,± 



> 2T . 



(84) 

Using the regularized Cooper pair function, we can eval- 
uate the Cooper pair size (Coherence length) 0, [Z^] : 



£{T) 



(0 + (r)|r 2 <Mr)) 
(Mr)\Mr)) 



(<z)|v 2 <Mg)> 



(0 + (q)|<M<?))A 



(85) 



FIG. EKbl) shows the Cooper pair wavefunctions 
4:TT 2 k F r(t)±(r) at T = T c in the BCS phase (G C /G R = 
—0.35). We can see that the wave function <t>+(r) oscil- 
lates with a period ~ 2?r/kp and decays exponentially 
c with the Pippard length 



k, 



1 



Ef ttAq 



(86) 



Ao ~ 1.764T C is the gap at zero temperature [44|. The 
Pippard length is the weak coupling (BCS) approxima- 
tion of the pair size at T = 0, i.e., £ p = ^ C (T = 0) 
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FIG. 5: The spectral density as a function of w + 2/j, (left panel), wavefunctions in the real space (middle), and the occupation 
numbers as a function of momentum k (right panel). From top to bottom, the attractive coupling increases as G c /Gr = —0.35 
(a), Gc/Gr = 0.5 (b), and G c /Gr = 1.35 (c); the system evolves as the BCS (a), the BEC (b), and the RBEC state (c). The 
inset of (al) is an enlargement near uj — 0, and the inset of (a2) shows the behaviour of the wavefunction at large length scale. 



in the BCS regime. The inset of FIG. EJlb) shows 
the magnitude of wavefunction actually decays exponen- 
tially at large scale r > £ c = £ C (T C ); note that the 
pair size is robust against the increase in temperature, 
£c£T c ) ~ £ C (T = 0) [7i|. Also it is consistent with 
|78| that the antiquark pair correlation is negligible in 



the BCS regime. How this situation changes as the 
coupling grows? FIG. E)j2b) shows the wavefunctions 
Air 2 kFr<f)±(r) in the BEC region. We first notice that 
the wavefunction 4>+{r) has no oscillating structure like 
that in the BCS. This is easily anticipated by the fact 
that the in-medium bound state forms in this regime and 
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that the bound state in the s-wave channel is unique; so 
the "wavefunction" should not have any nodes in the real 
space ■ This is in contrast to [43|, [3] where the pair 
structure at low density still has a oscillating structure; 
this fact shows that the in-medium bound state is not 
formed there. In this region, the Pippard length is not a 
good approximation of the pair size. In this regime, the 
pair size is roughly approximated by the s-wave scatter- 
ing length a s . This is because, neglecting matter effects, 
the bound state wavefunction in the CM. frame should 
behave as ~ e as r — ► oo. We will confirm these 

later. It is also notable that there is a slight deviation in 
the wavefunctions 4 > +{ r ) an d 4'-{ r ) m this regime reflect- 
ing the explicit breaking of charge conjugation by finite 
fx. In the RBEC phase, however, <f>+ and 0_ show almost 
the same behaviour because of \i c <C T c (see FIG.[Sj3b)). 

Finally, we discuss how the Fermi surface is destructed 
when the attraction is increased. To see this, we have 
plotted the occupation numbers in FIG. Hfcl) ~ (c3). 
From these figures, we can clearly see that the Fermi sur- 
face gets significantly broken once the in-medium bound 
state have appeared in the spectrum. Baryon number 
carrier in the BEC regime is almost the bosonic degrees 
of freedom as we can see also from FIG. [TJ(b). In the 
RBEC phase, the temperature is high T c ^> /i c , and thus 
fermions are again thermally active for the entropic rea- 
son we noted in the previous section. Also, the quark- 
antiquark asymmetry is small because of (fi c <C T c ). 

IV. THE EVOLUTION OF THE SOFT MODE 
DYNAMICS 

In this section, we construct the low energy effective 
theory for the pair fluctuation transport and study how 
the static and dynamic nature of the soft mode evolves 
as the system goes from the weak coupling BCS to the 



RBEC. In the context of a nonrelativistic BCS supercon- 
ductor, the effective theories for the fluctuating pair field 
both for (i) near T c and (ii) near T = were derived in 
[8C| where a diffusion type equation was found above T c . 
The dynamics of the spece-time variation of fluctuation 
in color superconductivity above T c was first investigated 
in pTi ] where the authors found a damped-oscillation type 
equation within the linear level. The effect of fluctuat- 
ing pair fields to the specific heat was also investigated 
in [8l| starting with the ansatz that the pair fluctuation 
is a propagating mode with a relativistic dispersion. By 
extending [47| along with the nonrelativistic calculations 
we will find the following three points: (i) The 
pair fluctuation in the BCS regime is a diffusive mode 
but with small propagating piece [13, [z^|. (ii) In the 
BEC regime, the fluctuating pair field is of a propagat- 
ing mode but with a nonrelativistic dispersion [76[. (iii) 
The fluctuation is also propagating in the RBEC region 
and the relativistic (an almost linear) dispersion is real- 
ized in a wide kinematical region. 

Let us start our analysis with first looking at the static 
part of the fluctuation. In the context of color supercon- 
ductivity, this was first done in [82| in the weak coupling 
BCS regime. The long wavelength expansion of the pair 
susceptibility near T = T c leads to 

[F^(0, P)} - 1 = + £P 2 + • • • , (87) 

where the mass-squared parameter ao and the stiffness 
parameter (or diffusion constant) c arc defined by 

_i ' (88) 
c = 4m^[r^(0,P)] | p=Q . 

The explicit calculation yields [83[ : 
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(89) 



We note that ao has no UV divergence, while c still has 
a logarithmic divergence in the last term in Eq. (|89")1 
which is sub-sub-dominant in the (/i/T c )-expansion; ac- 
cordingly, the integral is conversing in the weak coupling 
limit, while it is not in the strong coupling. We may 



simply regularize the divergence by sharp cutoff A. 

In the extremely weak coupling BCS regime where 
fi c 3> T c , we can derive some analytic approximations. 
First, we may safely ignore the contribution from anti- 
quarks, and can replace the momentum integral j j^prj 
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by N J_ x dE q with N 
at E, 



2jt- 



being the density of state Further explicit calculation leads 



q — fj,. We then arrive at the approximations: 
BCS _ 7C(3)fcf 



a 



BCS 




47V, 



4m 12ttVT c 2 



F N . 



Because r*f(0,P) 



P2 + i 



near T = T c . 



(90) 



ith 



T-T c 



being the reduced temperature, the static dis- 
turbance of the system restores at the "healing" length 

£giA7 1/2 , (91) 

where the Ginzburg-Landau (GL) coherence length £gl 
is defined by 



6 



GL 



(92) 



The healing length diverges as T — » T c + as a conse- 
quence of the second order phase transition. In the weak 
coupling BCS regime (/i ^> T c ), the GL coherence length 
is approximated by (substituting Eq. (|90|) into Eq. (|92|1 ). 



eBCS 



fcjr V21C(3) 



1271. 



0.739^, 



(93) 



a well-known relation between the Pippard length and 
the GL healing length in the nonrelativistic BCS theory 
We now look at the dynamic part of the problem. One 
must be careful in performing the low frequency expan- 
sion of the dynamic pair susceptibility near T c because 
there arises a spectral singularity in the kinematical re- 
gion oj/P < vp once the system goes below T c [8(|; vf is 
the Fermi velocity of quasi-quarks. This singularity is at- 
tributed to the local absorption (or emission) of the Gold- 
stone mode (phonon) by thermally active quasi-quarks. 
A simple expansion is allowed only for the region near 
the critical temperature T ~ T c where the gap is small 
lu A(T). In such region, we can naively perform the 
low frequency expansion [80j as 



[rJ$f(w,o)] = -du-d 2 u 2 + ■•• 

where d and d 2 arc defined by 

d = -^[r^(^o)]- 1 L =0 , 
d 2 = ^[r^>,o)rU. 



(94) 



(95) 



These quantities are complex in the BCS region (fi > m), 
while they become real in the BEC region (fi < m) be- 
cause the imaginary part of the dynamic pair suscepti- 
bility looks like 



im [rJ$f>,o)]- 
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(97) 



The weak coupling (/z 3> T c ) approximation again applies 
both to d and d 2 yielding the following analytic results. 



d BCS = iNo^- 
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(98) 



Note that d BCS is pure imaginary reflecting the particle- 
hole symmetry in the weak coupling limit [801 ] . This 
means that, to this order, the pair fluctuation obeys a 
diffusion-type equation in the BCS phase [76|, |8(j. In re- 
ality, there is a finite real part which is an order (T c / fi) 
suppressed from the imaginary part though. This real 
part arises from a particle-hole asymmetry and adds a 
small propagating nature to the soft mode; as a con- 
sequence, the fluctuation in the BCS region becomes 
a damped- oscillation mode (47| or an overdamped mode 
[76| . The coefficient of second-time derivative, d 2 , is also 
non- vanishing. This is another source of a propagating 
nature of the soft mode. 

So far, we have concentrated on the kinematical as- 
pect of the pair fluctuation. It is sometimes important to 
know about the interaction between soft modes to inves- 
tigate the transport properties of system. For example, 
the shear viscosity is proportional to the mean free path 
(^mfp) of soft modes in the classical level. This requires 
an information about the collision between soft modes. 
We here look at the interaction between soft modes. 

The partition function is factorized as Zoe~ Sc!f where 
Zo is the partition function for free quark gas, and S e g is 
the effective action for fluctuation. We need to go beyond 
the gaussian approximation for S e g. Up to the quartic 
order in A, we have 
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S cff [A,A*] = TJ2J 7S3 [^(Mn,P)] "lA^iQ^P) 
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(99) 



Here numbers {1, 2, 3} arc labeling sets of frequency and 
momentum {(ifijvij -Pi), (i^N2, P2), (i^m, P3)}, and 
Si, 2,3 is a shorthand of T 3 J2ni,n 2 ,N 3 I (2T+ 7 (TTp 
To derive the low energy effective action, we set 
A J? (if2jv) = for iV = ±1,±2, and only con- 
sider the fluctuation mode with the lowest Matsubara 
frequency. Thus A has no dependence on r. Then 
the quadratic term of S c ff is casted into the form 
P J dxf^g 1 (A(x), A* (x)) with the kinetic part of the local 
effective (classical) free energy density defined by 



A n (x) 



= E„ MA ?) 0r)| 2 + ^\V x A v (x)\ 2 ] 



ax 



a o T T Tc is the mass parameter and c is the diffusion 



constant introduced before. The quartic term also can be 
written in the form j3 J dxf$(A,A*) with 
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(100) where b = 60,0,0 = 



0,0,0 



is calculated as 



(101) 
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(102) 



To derive this, we have taken care the lowest term omit- 
ting the gradient terms. This corresponds to the soft 
limit of the two body scattering where all the incoming 
and outgoing momenta vanish. The weak coupling BCS 
(fx ^> T c ) approximation again applies to bo resulting in 
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(103) 



Collecting all the above results, we arrive at the 
following expression for the local (but coarse-grained) 
Ginzburg-Landau (GL) free energy functional (F c g = 
J dx(f^ + )) which describes the long wavelength 
(static) excitation of soft mode: 



-Foff [A, A*] = Jdx £ Ja T |A r) (aO| 2 + ^| V^A^)! 2 ] + ^^\A n (x)\ 2 \A 6 (x) 



(104) 



On the other hand, we have already derived the time dependence of the excitation at long time scale. Combining the 
nonlinear term in A(t, x), we obtain the following dynamic equation for the fluctuation transport: 



-iddt 



d 2 d 2 + aT -^V 2 +&oEl A ^( ^ > ; 
l 



A v (t,x) = 0. 



(105) 



It can be written in the compact and familiar form 

<5F eff [A,A*l 



As for extremely low frequency modes with (u> <C 



[-iddt +d 2 d?] A v {t,x) 



SA*(t,x) 



(106) 



17 



> 9 9 
-9 



6 nnn 2™. X > 

0,0,0 _ s6_ s ^ 

2d 2 M fl _ ¥ X 



-Jd\" 



■JdA"' 



FIG. 6: The microscopic process which contributes the boson- 
boson repulsive force in the BEC region. The thin line ex- 
presses propagation of fermion, while the external dashed line 
represents incoming/outgoing boson. \fd is the wavefunction 
renormalization. The arrow indicated on the line represents 
the flow of baryon charge. 



min. (T c , \m— fx\)), we may simply ignore the second term 
(including di) in the left hand side of the above equation. 

In the BCS regime, d is imaginary dominant and is 
pure imaginary in the weak coupling limit as shown in 
Eq. (|98j) . In this regime, Eq. (|106[) describes how the 
fluctuation about the thermal equilibrium relaxes; the 
thermodynamic restoring force to the equilibrium is then 
given by -^^ ^l^ 1 ■ By rescaling the field by 

ty(t, x) = y / cA(t, x), we find the standard expression for 
the time-dependent Ginzburg-Landau (TDGL) equation 
in the BCS regime: 
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(107) 

The pai r fluctuation is diffusive as in the nonrelativistic 
case j76l l80l ] ; ignoring the nonlinear term, the mode with 
wavenumber k decays exponentially 
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with rfc being the relaxation time 
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(108) 



(109) 



GL 



The relaxation time diverges as k — > and T — > T c + 
reflecting the critical slowing down. The TDGL equation 
for fluctuating diquarks is derived in the linear level [8^ | . 

Now we discuss the fluctuation transport property in 
the BEC regime. In the BEC regime where fi < m is 
realized, d in turn becomes pure real as shown in Eq. (|96[) . 
Again by ignoring di term looking at the dynamics of 
sufficiently low energy excitation, and by rescaling the 
field as ^(t,x) = y/dA(t,x), we arrive at 
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It is obvious that this is the Gross-Pitaevskii (GP) equa- 
tion [84, 85, 86] which describes the low energy excitation 
above the Bosc-Einstein condensate after identifying 



Mb 



"1 ' Mb 



„ d 
2m-, 

c 



M h 



where ^h, Mb, and a s ;, are the effective chemical poten- 
tial, the effective boson mass, and the effective scattering 
length for s-wave boson-boson scattering. In Fig. [6l we 
show the schematic Feynman graph which contributes 
to a s b- When the temperature goes below T c , Hb be- 
comes positive indicating the instability to formation of 
the Bose-Einstain condensate. 

We see that the fluctuation becomes propagating in 
the BEC regime as in the nonrelativistic fermion system 
[76| . This is because a gap (2m — 2/i) appears between 
the continuum excitation and u> = due to the formation 
of bound state. Therefore, the low energy fluctuation 
about the equilibrium cannot decay in this region; the 
fluctuation with wavenumber k can decay only via the 
nonlinear term, i.e., the scattering among soft modes. 

We now study the situation \du>\ c^w 2 ; we will see 
that, in the RBEC phase, there actually exists such an 
"intermediate" frequency regime d/d 2 ^ \w\ where the 
expansion is still valid (cu -C cu c = min. (T c , \fx — m\)). In 
such regime, we find di becomes pure real and positive. 
By rescaling field ^(tjX) = y/d~2A(t,x), we obtain the 
relativistic Gross-Pitaevskii (RGP) equation (or simply 
the Klein-Gordon equation with $ 4 interaction) (87|: 



d?-v^l+AP + \J2\*dt^)\' 



%n(t,x) = 0, 



(112) 

where v s , M 2 are the velocity and the effective mass of 
soft mode. A represents the two body repulsive force 
between soft modes. They are 
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(113) 



When the temperature goes below T c , M 2 becomes nega- 
tive. The soft mode then becomes tachyonic without non- 
linear terms, which signals the instability to the forma- 
tion of nonzero Bosc-Einstein condensate (^(0, 0)) ^ 0. 

To see which region of (u>, k) -space is governed by the 
RGP (or GP) equation, we look at the dispersion relation 
of the fluctuation in the (R)BEC region. The dispersion 
is given by the solution of 
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At the critical temperature (ot = 0), we find 
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FIG. 7: (a) The static GL parameters as a function of G^ 1 . ao, c, and 6 are the mass-squared parameter, the diffusion 
constant, and the interaction parameter. Oq CS , c BCS , and 6 s are the weak coupling analytic approximations (/i c 3> T c ) of 
ao, c, and 6. (b) The dynamic GL parameters, d (efe) is the coefficient for the first (second) time-derivative of the TDGL 
equation. lo c = min.(T c , |/x c — m\) is the limit energy such that the low energy GL expansion is valid for w C w c . 



For a sufficiently small wavenumber k -C 



2M, ' 



(116) 



where Mb is defined in Eq. ()111[) . This is the dispersion 
of nonrelativistic particle with velocity v = k/Mb- In the 
opposite case, the dispersion is approximated as 



(117) 



with v s defined by Eq. (|113[) . This is the dispersion of 
phonon with velocity v s . We note that as di becomes 
large such that cfo 3> cPm/c, the momentum region where 
the dynamics of fluctuation is described by the RGP 
equation becomes wide. This is the case when we in- 
crease the attraction between quarks as we will see later. 

In FIG. [UJa) and (b) show the numerical results for 
the static and dynamic GL parameters from the weak 
coupling BCS regime to the ultra-strong coupling RBEC 
regime. The weak coupling approximation for these pa- 
rameters are also indicated by several thin dashed lines. 
One can see the drastic change going from the BCS to 
the RBEC particularly in the dynamic GL parameters. 
From these results, we shall discuss the character change 
of the soft mode, and transport properties. 



A. How does the static property of the soft mode 
evolve with the crossover? 

We first focus on the GL coherence length Eq. ([5^)1 . 
As noted, the static disturbance to the system restores 
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FIG. 8: The GL coherence length £gl and the Cooper pair 
size £ c evaluated on (/j,c(Gr),T c (Gr)). The Pippard length is 
shown by the line indicated by a a is the scattering length 
for s-wave quark-quark scattering. These are plotted in the 



unit of the mean inter-fermion distance N t 



-1/3 



at this healing length. Also the anomalous heat capacity 
(5C a per unit volume is controlled by this length as 
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(118) 



Furthermore, this quantity is known to appear in the 
intensity of the diverging conductivity at T = T c (88l . l89l . 
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B. What about the dynamic nature of fluctuation? 



6 



'1/2 



GL 



(119) 



where g is the gauge coupling. 

In FIG. [3 we show the GL coherence length £ql as a 
function of G^ . The Cooper pair size £, C (T C ) and the 
Pippard length £ p are also shown. In the weak cou- 
pling region, all these scales exhibit almost the same 
behaviour. In particular, the weak coupling "universal 
relation" £gl/£ p = 0.74 (see Eq. (|9"3"j) ) is satisfied in 
good accuracy. However, once one goes into the BEC 
region, these scales start to deviate. First, the pair size 
£ c becomes smaller than the mean inter-fermion distance, 
which supports the picture of independent bosons. In this 
regime, the pair size £ c is approximated by the scatter- 
ing length a s rather than £ p as noted before. In contrast, 
the GL coherence length £ql takes a minimum between 
the BCS and BEC regimes as nonrelativistic calculation 
[79l | . On the other hand, the behaviour of £gl is ap- 
proximated by the Pippard length £ p ~ vp/Ao which is 
defined in the weak coupling regime though. This fact 
suggests that, the gap parameter not only has the role 
of the order parameter but also plays a physical role to 
determine the response to the static disturbance. 

Ginzburg-Levanyuk region: The above derived low 
energy effective theories may fail in describing the dy- 
namics of soft modes when T is very close to T c , i.e., 
they cannot be applied in the critical region. In this crit- 
ical region the thermodynamic quantities diverge with 
the anomalous power law. We can estimate this region 
by Ginzburg-Levanyuk criterion [9l], [92j. We first de- 
fine the order parameter averaged in the sphere Vt = 
47t(£gl/\/^\) 3 ! a characteristic size of static fluctuation: 



Vt 
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Then we estimate the magnitude of the thermal average 
of fluctuation in the gaussian approximation. 
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If aT(^'* ) x I , ) ) )th becomes smaller than the quartic term 
||^| 2 |^| 2 , then the mean field approximation fails due 
to higher terms in the soft mode action; this condition is 
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When the temperature is in the above range, the 
Ginzburg-Landau approximation fails to describe the 
anomalous divergence of the thermodynamic quantities 
like the heat capacity (47[. We must note that the criti- 
cal region estimated by above formula is not small in the 
RBEC regime, which means we need to go beyond the 
quartic approximation in order to describe the soft mode 
there in the quantitative level. 



We now look at the dynamic property. We first ignore 
the effect of di looking at the low frequency regime. Then 
the complex pole (w po ic) of the dynamic susceptibility in 
the TDGL approximation should be given by 
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where the effective (complex) mass M e g, the relaxation 
time t r , and the oscillation time t° are defined by 
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In the BEC region, M e s is real and coincides with effec- 
tive boson mass Mb defined in Eq. (|111[) . For the stability 
that the fluctuation about the thermal equilibrium is not 
tachyonic, Imcj po i c < should be satisfied. This requires 
both T > T c and Ime? > 0, and our numerical result 
shows that the latter is actually satisfied (see FIG.[7Jb)). 
The fluctuation with a wavenumber k behaves as 
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the fluctuation looks like A^(i, fc)|r=T c = 
3 + const.) where if = 2 Ji^li I -j 7 _ and 

-p — R C I/ ~ arc t nc diffusion and propagating 
times. Thus the complex effective mass determines 
whether the fluctuation with a finite momentum k is dif- 
fusive or propagating. In contrast, t r and r° determine 
how the long wavelength limit of the fluctuation behaves 
for T T r . In this case, the fluctuation looks like as 



(120) A v(t,0) 



sm(t/T°t r 1 + const.) where t r and 



t are the relaxation and oscillation times. 

FIG. [9]Ja) shows the complex effective mass M e g as 
a function of Gr, and (b) shows the relaxation time 
t R (Gr) and the oscillation time t°(Gr). We can see 
from figure (a) that the fluctuation is diffusive (t r ^> ) 

(121) in the BCS regime, while in the BEC regime, the pair 
fluctuation is of pro pagating (t r <C t®) because of the 
reality of boso n 17611 . In contrast to the nonrelativistic 
calculation [761 177|. |M e ff| ^ 2m in the strong coupling 
limit. It decreases due to the binding effect as the cou- 
pling becomes stronger. Also the long wavelength limit 
of the pair excitation is of overdamped in the BCS side 

(122) (t r t°) while it is oscillating mode in the BEC side 
(t r = 00) because of Imd = 0. This is again attributed 
to the bound state gap. 



C. Relativistic Gross-Pitaevskii equation 

We have shown that there may exist the three distinct- 
(u>, k) regions where the fluctuation is described by the 
TDGL, GP, and RGP equations. FIG. flTJt a) shows where 
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FIG. 9: (a) The complex effective mass as a function of Gr, which determines the time dependence of the fluctuation with 
finite k at T = T c . (b) The oscillation time r° and the relaxation time r R for the fluctuation in the long wavelength limit. 



1.0 




-0-5 0.0 0.5 1.0 1.5 '-0.5 0.0 0.5 1.0 1.5 

GJG R GJG R 

FIG. 10: (a) The three regions where the fluctuation is described by TDGL, GP, and RGP equations. The low energy expansion 
is allowed for the shaded area where ui < uj c = min.(T c , \fi c — m\). (b) The group velocity of the fluctuation (boson), duk/dk, 
evaluated along the light dashed arrow in (a). 



in the w-rcgion these three equations govern the dynam- 
ics of soft mode. The shaded area satisfies the condition 
cu < cu c = min.(T c , |/i c — m\) where the low energy ex- 
pansion is valid. As in the nonrelativistic case [7g|, the 
validity breaks down in the very vicinity of the BCS-BEC 
crossover due to |/i c — m\ — > 0. The solid line indicated 
by wgp/rgp divides the w-space into two pieces, one in 
which the GP equation applies, and the other where the 
RGP equation applies. We determined this boundary by 

WGP/RGP = TT, . ( 126 ) 



which is evaluated by equating ^jk 2 in Eq. (|115| to 
unity; if u) 3> wgp /rgp the dispersion of the fluctuation is 
approximated by Eq. (|117p . and in the opposite case with 
lu <C wgp/rgp, it is approximated by the nonrelativistic 
dispersion, Eq. (|116p . We see that the w-region where 
the RGP equation describes the fluctuation rapidly grows 
when the system goes into the RBEC region. This is 
obviously because boson is light in the RBEC region. 

FIG. fTOl b) shows the group velocity of the fluctuation 
mode, i.e., v s — (with Uk defined by Eq. (| 1 1 5|) ) eval- 
uated along the dashed arrow in figure (a). The velocity 
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FIG. 11: The fermion-fermion scattering length a s = 
tuGr/Ait, and the boson-boson scattering length a s b- The 
arrow indicates 2, the strong coupling limit which the ratio 
a s b/a s approaches in the nonrelativistic calculation [7r3 ]. 



of fluctuation is smaller than the Fermi velocity vf and is 
almost independent of Gr in the BEC region. It rapidly 
grows at G c /Gr ~ 1 where the system goes into the 
RBEC phase; the speed v s becomes comparable to the 
speed of light; this is also consistent with the fact that 
the boson is light in this regime. 



D. Interaction between soft modes 

As we have discussed, there are multi-body interaction 
among soft modes. To the lowest order beyond the gaus- 
sian approximation, we found the repulsive two-body in- 
teraction which causes a binary collision. From Eq. 
the scattering length for this collision is given by 



a>sb 



27rcrf' 



(127) 



which determines the scattering cross section in the soft 
limit (s — > 4M^ ) in the CM. frame as 



0-6 



,2 1,2 



2 d 2 ' 



(128) 



In Fig. 111! we show the fermion-fermion scattering 
length ma s and the boson-boson scattering length ma s b 
as a function of G^ 1 . Also its ratio a s h/a s is depicted. 
When G^ 1 is increased beyond the BCS/BEC bound- 
ary, the boson scattering length first decreases, i.e., the 
weaker boson-boson repulsion in the stronger fermion- 
fermion attraction uS . This is because the boson-boson 



interaction is microscopically induced only via the dis- 
sociation of boson in the intermediate state as shown 
in Fig. [SJ it becomes harder to dissociate the boson in 
the stronger coupling. When G^ 1 is further increased 
and the RBEC is approached, the boson-boson scatter- 
ing length rapidly increases, and the nonrelativistic limit 
value of a s b/a s , 2 [7^|, is significantly exceeded 1 . This is 
because the temperature is large (T c > 2/x c ) in the RBEC 
regime, which favors dissociation. 



E. Smooth change of transport properties 

We here discuss some of the transport properties on 
the basis of the dynamic equations obtained in the pre- 
vious section. In the nonrelativistic superconductor, it 
is known that the fluctuati ng p air fields affect the elec- 
tric conductivity above T c [88l . l89l . l90j . In the context 
of two-flavor color superconductivity, it is noted in [45[ 
that these soft modes may affect not only transport co- 
efficients but also the dilepton spectrum. 

We here focus on the shear viscosity near the criti- 
cal temperature and derive approximate formulas for the 
viscosities in three regions. To do this, we employ the 
semi-classical Kubo's formula and also the Boltzmann 
equation for soft modes in the sprit that the short range 
quantum effects are already integrated out in the low- 
energy coefficients in the soft mode action. For this to 
be correct, soft modes should be only dilutely distributed 
in the phase space such that they can be described by the 
classical effective theory. Technically, this means that we 
only take the "thermal average" of the soft mode dis- 
tribution in the Kubo's formula. We will make a short 
statement regarding the validity of this treatment in the 
end of this section. 

Shear viscosity due to the soft mode diffusion in 
the BCS regime: We now look at the precursory soft 
mode contribution to the shear viscosity, r] s , in the BCS 
regime. In the rotationally symmetric system considering 
here, the shear viscosity is formally given by the quantum 
Kubo's formula (9f|, 

n = lim - fdx f dte i( - u+i5 ^([f xv (x,t),f xv (0,0)}), 
u J Jo 

(129) 

where denotes the traceless part of the energy- 
momentum tensor. First thing we have to do is to derive 



1 In contrast to the approach taken in [7611 . the study of 4-body 
Schrodinger equation produces as ~ 0.6a s near the unitarity 
[93ll . In addition, the renormalization group approach yields the 
similar value also in the strong coupling regime |94| . Thus it is 
clear that, for the quantitative understanding of the dimer inter- 
actions, we need to go beyond the current approximation taking 
into account how higher energy 2-boson processes renormalize 
the low energy scattering length. We defer this to future study. 
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the classical approximation (u> <C T) of the Kubo's for- 
mula Eq. (|129p . From the spectral representation of the 
real part of Eq. ([1291) , we find for u <C T, 

Re 77 = p J dx J dt(f xy (t,x)f xy (0,0)}. (130) 

Note that we have no longer the commutation in the 
bracket. We now replace T xy by the classical stress tensor 
and the expectation by the thermal average looking only 
at the contribution from soft modes (with u> <C T) to the 
viscosity rj s . For this purpose, we first extract the soft 



mode contribution to the stress-momentum tensor which 
we denote by STij . The stress arising from soft modes is 

STij{t,x) = -%/ off (A,A*) 

+^[d i A* rl (t,x)d j A v (t,x) + (i^j)] 

(131) 

This is followed by calculation of the free energy shift 
SF g = — J dx8Tij-J=^-Ui{x) under the small distortion 

i(x). The shear viscosity from fluctuat- 
ing pair fields is then given by 



/Him f°°dt^- fvAVA*e-^ A ' A "^T xy (t,p)5T xy (0,-p) 



(132) 



Here (■ • -)th denotes the thermal average at t = 0, i.e., 
-^3— J T> ADA* (■■■). The diagrammatic interpretation 
for the above formula is given in Fig. 1121 To proceed 
further, we use the standard ansatz that the time evo- 
lution of the fluctuation is determined by the classical 
TDGL equation in the gaussian approximation; 



U qi tA q Q 



fa(q 2 +€GL*-)A 



(133) 



The gaussian approximation to F B s[A, A*] is also es- 
sential in order for the replacement of the expectation 
of A 4 by the product of the two body Green's func- 
tions in Eq. p32|) to be valid. By taking the identity 

(A^ A I? * )th = ~~~~ o r-?=3 . into account, we obtain 
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where we have introduced the UV cutoff q c 



= k F ). 

(134) 
and set 



kp. Unlike the conductivity which diverges as 



|89| , the shear viscosity is regular. 



T^T C like a ' 

Also the shear viscosity is proportional to the soft mode 
relaxation time; the transport coefficients are always pro- 
portional to the relaxation time (mean free time) of the 
quasi-particle according to the relaxation time approxi- 
mation to the Boltzmann equation. Because of the diver- 
gence of the relaxation time t r — > oo at the BCS-BEC 
crossover point [i = to, this formula fails to describe the 
shear viscosity in the BEC regime. We thus try the other 
approach to derive the shear viscosity in the BEC region. 

Shear viscosity via boson binary collision in the 
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FIG. 12: The microscopic process which contributes to the 
shear viscosity in the BCS region (/x c > m). The same graph 
expresses the Aslamazov-Larkin type correction to the con- 
ductivity if the vertex Tij(t,x) is replaced by the paramag- 
netic current [47]]. In our TDGL approach, a phonon couples 
to two soft modes by the TDGL coefficient (c/4m). 



BEC regime: In the BEC regime, the soft fluctuation 
about the equilibrium cannot decay via the quantum dif- 
fusion because of the bound state gap. In reality, how- 
ever, the soft mode with momentum k decays via the 
nonlinear term in the GP equation, i.e., the binary col- 
lision of soft modes. Here, we simply apply the result 
from the classical Boltzmann equation for dilute gas of a 
particle with a hard sphere r = y/2a s b- [96;, 9~7l, l9a|. 



5tt 



n so ft \M e ffV)thlmfp 



— 0.35 X n S oft(M e ffu)thimfp. 



(135) 



n soft is the soft mode density, and i m f p is the mean free 
path. (M e {fv)th is the average momentum of thermally 
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excited soft modes. These are given by 

^mfp 



n ao it<J s b n soit m 2 bl'' 



(M cS v) th = X 7 



M r ffT / mdT 



(136) 



2ir V C7T 
Combining the above results, we arrive at the formula 

5^2 1 5ir 3/2 c 3 / 2 d^ 2 ^T 



256 X T a 2 b 32J2 m 3 / 2 b 2 ' 



(137) 



This equation relates the low energy coefficients in the 
GP equation to the shear viscosity. We note that this 
formula is valid only for the classical (dilute) limit, i.e., 

"softa^, n soft -C 1, where At = J m^t 1S tne thermal 
dc-Broglie length, the quantum radius at T. 

Shear viscosity via binary collision between light 
bosons in the RBEC phase: The nonrelativistic for- 
mula we have derived above fails in the RBEC regime 
because of the lightness of boson. The qualitative be- 
haviour of the shear viscosity is 0, Il00j | 



Vs ~ l m fpV(e + P)th ~ lmfpVT{s) t h, 



(138) 



where s denotes the entropy density, and v is the average 
velocity of soft modes. In the above formula, we ignored 



the chemical potential dependent part of the pressure, 
i.e., P = fin + Ts — e — ► Ts — e with e being the en- 
ergy densit y. T he relativistic thermal de-Broglie length 
is given by jlOl| 



T 2/3 
~f~' 



(139) 



and therefore, the soft mode dynamics is described by 
the RGP equation we derived in the previous section, v 
and the scattering cross section for a collision between 
soft modes is estimated as 



4md2 



<7 sb (s ~ T 2 ) 



A 2 



T 2 T 2 d\ ' 



(140) 



with A being the A 4 coupling (see Eq. (|113p ). We obtain 



v Ts 

<7sb Tlsoft 



(141) 



where we adopted the simplified notation s = (s)th for 
the equilibrium entropy density. n so ft again denotes the 
density of thermal soft modes. We further simplify the 
equation by estimating s/n so f t in the present gaussian 
approximation for the soft mode action. The thermody- 
namic potential coming from soft modes is 



SSI = M^R T J2 N J $p Z„=± log (-d 2 (iSl N + an) 2 + ^ + a T ) = SSl + + 5(1-, 
where SSl± represents the bosonic (antibosonic) contribution to the partition function: 



(142) 
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By differentiating the difference (8SI+ — SSI-) with respect 
to fi, we find the following formula for the soft mode 
density in the limit T — > T c and T c 3> /it, 



_ c . (d 2 m\3/ 2 16C(3) N C (N C -1) T 3 

«S0ft — { c ) ^2 2 1 C ■ 



(144) 



Also, we take the same limit in Eq. (|142p to find the 
Stcphan-Boltzmann pressure 



SSI = 



_ N C (N C ~1) 



tin 



(145) 
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Then by differentiating this with respect to T c , we find 



, v _ N C (N C -1) Air 2 
\ s hh - 2 45" 



(146) 



Combining all the above, we obtain the following para- 
metric dependence of the shear viscosity 

Vs ~ \/d^v—, (147) 

CTsb 



Then, adjusting the prefactor to the result of the rela- 
tivistic Boltzmann equation [98[ , we end up with 



Vs = 



3 / c T c 
10-7T V 4m a s h 



3 c 1/2 4 3 
207rm 1 /26 2 , c ' 



(148) 



Summarizing our results, shear viscosities arising from 
soft modes in the three regimes arc 



jfc(BCS) 
Vs(BEC) 
77 S (RBEC) 
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FIG. 13: (a) The shear viscosity to entropy ratios from soft modes, r/f cs , r/f EC , and 77^ BEC as a function of four-fermion 
coupling G^ 1 (thin solid lines). The dashed lines indicated by ?7/ CS 77/ EC , an d 77^ BEC correspond to the shear viscosities from 
fermions in BCS/BEC and RBEC regimes. The total viscosity, the sum of the fermion- and boson- originated viscosities, are 
shown by bold solid line. The "KSS bound", 1/4-7T, is shown by the horizontal dashed line. At the unitary limit, the cross 
section for the low energy s-wave fermion scattering diverges as ~ 1/k 2 with k being the relative momentum, (b) The shear 
viscosity to quark number ratios. The dashed line indicated by "GSZ bound" is the minimum bound, l/6n, proposed in [68( |. 



The system has fermion's degrees of freedom as well as 
soft modes. As for the viscosity from fermion-fermion 
binary collision, we estimate that in the nonrelativistic 

if^T^T with * bcing the avera S e 
momentum of fermions, and that in the RBEC regime 

by Vf 



region by rjf 



- . We also evaluate k simply by k 



40tt 2 m?o?. 



max. (/1 2 — m 2 , ^^), i.e., we take the larger of (a) the 
matter momentum due to the Pauli-blocking, and (b) the 
thermal momentum. 

In FIG. H31 we have examined the above-derived for- 
mulas for shear viscosities in three regimes. FIG. Hffl a) 
shows the shear viscosity to entropy density ratio as a 
function of Gn . This quantity is related to the sound 



attenuation length by T s 



=^ with C being the bulk 



viscosity. For a hard sphere particle gas, ( -C rj both 
for nonrelativistic T -C m and relativistic T 3> m situ- 
ations [9^, [13, [HI • For the fluid-dynamic picture to be 
valid, this length must be much smaller than the time 
scale of the fluid-dynamic evolution of system. From the 
figure, we can see that in the weak coupling BCS regime, 
the total viscosity is dominated by contribution from the 
fermion binary scattering, and that from the soft mode 
diffusion is negligible 2 . As the coupling G^, 1 becomes 



2 We here estimated the total viscosity just by adding up the 
fcrmionic viscosity and the bosonic viscosity. In the kinetic the- 



larger, ?t^ cs monotonically decreases with the behaviour 
~ a 2 ~ 1/\Gr\ 2 ; because the cross section for the fermion 
binary scattering monotonically increases and this effect 
on the shear viscosity overcomes the effect of the de- 
crease in the average momentum k. In contrast, r]f 
monotonically increases because the relaxation time of 
the soft mode, trT Ci becomes large. As the unitary limit 
G^ 1 = 0~ (a s = —00) is approached, r]J cs rapidly de- 
creases and crosses the "KSS bound" (l/47r), the mini- 
mum bound of rj/s conjectured out by Kovtun, Son and 



ory, the viscosity of the gas mixture takes much more compli- 
cated form even in the first approximation of the Chapman- 
Enskog equation where the deviation of the distribution func- 
tion is expanded up to the first order in the Sonine's polyno- 
mials. Let us here quote the result of [97j (see Eq. 12.5.1); 
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with ni (JI2) being the number density of each particle in the 
gas mixture, fix (fJ.2) being the viscosity of the individual gas, 
and mi {1^2) being the mass of each species. A and E depend 
on the microscopic detail of interaction between cross species, 
and in most cases, E is inversely proportional to <r 2 2 where 1T12 
is the collisional cross section between cross species. Then, if 
there is no interaction between cross species, E is 00 and the 
above formula for r? m j x reduces to the simple sum of each vis- 
cosity, i.e., liniE^oo ?j m i x = 771 + r?2- In our case, it is simply 
assumed that there is no interaction between boson and fermion 
as a first step. Of course, further investigations are needed to 
judge if this treatment is quantitatively appropriate. 
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Starinets [ll| . However owing to the viscosity due to soft 
modes, the sum of the viscosity ratios, r] s + rjf, does not 
lower the bound. 

When the crossover point G c /Gr ~ 0.068 where [i c = 
m is closely approached, ?7^ cs rapidly diverges due to 
the divergent relaxation time (tr = oo) at /i c = m. As 
noted, however, this is not physical for the following two 
reasons: (i) The soft mode with k actually relaxes via the 
non-linear interaction caused in the TDGL action, (ii) 
The low energy expansion fails in the vicinity of [i c = m 
because of the absence of expansion scales. In short, the 
singularity in the crossover point is attributed both to 
our gaussian approximation to the soft mode relaxation, 
and to the breakdown of the low energy expansion [76| . 

In summary, we found the following facts. (1) Soft 
mode may play important role to the shear viscosity near 
the unitary regime. In fact, the model calculation per- 
formed here shows that it makes the minimum bound 
for shear viscosity proposed in remain valid. On 
the other hand, it is completely dominated by fermionic 
processes in the weak coupling BCS regime. (2) Shear 
viscosity takes almost minimum near the unitary limit 
in the present model calculation. It is close to the mini- 
mum bounds [ll|,[68j|. It is also noteworthy that the bulk 
viscosity vanishes in the unitary limit [70j. These sug- 
gest that the intermediate of the BCS/BEC crossover is 
rather liquid- like as conjectured in [68| . 

Finally we must note that our ansatz that the short 
range quantum correlation is completely incorporated 
in the coefficients in the low energy soft mode action 
may not be good approximation in the BEC and RBEC 
regimes. In fact, the ratio of the thermal de-Broglic 
length of soft mode to the average distance between soft 
modes is not smaller than the unity, taking around 2^4. 
Further theoretical investigations are indeed needed in fu- 
ture to clarify how secondary quantum corrections affect 
the shear viscosity. 



V. THE DENSITY DEPENDENCE AND THE 
QUANTUM FLUCTUATION 

We here again look at some static aspects of crossover. 
In Sec. IV Al we ask how the characteristics of the criti- 
cal temperature as a function of G^ 1 depends on quark 
number density which we have fixed so far. We then dis- 
cuss whether quark matter with the standard choice of 
diquark coupling exhibits the crossover. In Sec. IV Bl we 
look at how large the effect of quantum fluctuation in the 
"boson sea" , which we have neglected so far. 



A. The density dependence; what about the 
ultra-relativistic limit of the crossover? 

So far we have worked with a fixed Fermi momentum 
k F = 0.2m. We now study how the BCS/BEC/RBEC 
crossover will be affected by the change of quark den- 
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FIG. 14: The /^-dependence of the critical temperature as 
a function of G c /Gr. From bottom to top, the Fermi mo- 
mentum increases as kF = 0.2m, 0.4m, •• ■ ,2.4m. The large 
point located on each critical line represents the BCS/BEC 
crossover point where fj, c = m holds, while the large cross cor- 
responds to the BEC/RBEC boundary which we estimated 
by T c = 2/i c . The large square put on the horizontal axis 
indicates the standard choice of the diquark attraction, i.e., 
§ ~ |G S with G a = 2.17/A 2 set to reproduce the dynamical 
quark mass M q = 400 MeV in vacuum. 



sity. In FIG. [14j we show the fci?-dcpcndcncc of T c as a 
function of the four- Fermi coupling, G^ 1 . As the system 
becomes denser, the critical temperature increases mono- 
tonically 3 . At the same time, the BCS/BEC crossover 
boundary (jx c = m) indicated by the large point shifts to 
a larger coupling. As noted before, this is because the 
Pauli-blocking prevents the formation of bound boson in 
medium. In addition, the characteristic change in the 
shape of (T c , G^ 1 )-relation gets somewhat smeared and 
seems to vanish for k F *i m. Moreover, a typical diquark 
coupling strength, | of the scaler coupling G S A 2 = 2.17 
chosen so that the dynamical quark mass 400 MeV is re- 
produced at vacuum, is located below the unitary cou- 
pling G^ 1 = (see the large square on the horizontal 
axis). From these points, we conclude that quark mat- 
ter only with light flavors hardly exhibits the BCS-BEC 
crossover. We need a somewhat exotic situation where 
(i) the in-medium constituent quark mass is compara- 
ble to [i and (ii) the (qq) attraction between constituent 
quarks is stronger than that in the (qq) channel in or- 
der for the BCS-BEC crossover to be realized in possible 



3 One may think this statement is incorrect because the critical 
temperature seems to be shifted downward in the RBEC regime. 
This is not true; the density dependent normalization, Ep, is 
making T c at higher kp look smaller than those at lower densities. 
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quark matter inside compact stars. 



B. How does the quantum fluctuation affect the 
crossover? 

We here examine how large the quantum fluctuation we 
have ignored so far can be. As we discussed in Sec. Ill Bl 
the gaussian fluctuation consists of two parts, i.e., the 
Nozieres-Schmit-Rink correction 
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and the "quantum fluctuation" which remains finite as 
T — > in the presence of finite fj, 
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w-integral in O q fi is finitem, while the remaining momen- 
tum integral is quadratically divergent and thus we need 
to introduce a new three momentum cutoff As ■ In princi- 
ple, there exists no a priori relation between the fermionic 
cutoff A and the bosonic cutoff A B - But it is reasonable 
that Ab < A because the short range quantum effects are 
already taken in part to the bosonic effective lagrangian. 
Here we shall try two choices, Ab = 0.2 A and As = 0.3A. 

Before going into discussion of numerical results, let us 
little closely look at the structure of quantum fluctuation 
to see its physical meaning. In order to grab an intuition 
into this term, we try to estimate the quantum fluctua- 
tion by the bound state approximation for the in-medium 
phase shift; that is 



rd(uj < -Ef p - 2/0, 
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where E^'p and Eg'p are the boson and antiboson energy 
dispersions. Then substituting this into Eq. (|151|) yields 
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where the three momentum cutoff Ab is introduced. Now 
it is clear that this correction is due to the vacuum fluc- 
tuation (quantum Casimir pressure for boson). It differs 
from fermionic one by minus sign. For example, the chi- 
ral condensation energy at T = can be written by the 
integral of the zero-point energy shift as 



n xSB = -2N c Nftj^ 
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We conclude that the quantum fluctuation is interpreted 
as the "vacuum" energy coming from the in- medium shift 



of the boson (antiboson) dispersion. In the BCS regime, 
the above argument will be slightly modified because 
there is no stable boson in the spectrum. However, even 
in the BCS regime, it can be understood in the same way; 
it is the vacuum fluctuation coming from the in-medium 
spectral shift. 

Differentiating Eq. (|153[) with respect to fi results in 
the following expression for the number density from the 
quantum fluctuation. 
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Let us recall here that in the bound state approxima- 
tion, the number density coming from the ordinary ther- 
mal fluctuation is given by Eq. (|66p where we see that 
the "effective charge" of (anti)boson deviates from (±)2. 
From these points, we can see that a small fraction 
^^B(B)p/ °f c l uar k number charge in the thermal fluc- 
tuation escapes to the vacuum sector. In the nonrcla- 
tivistic situation where m 3> (|^ — m\, T c ), we can ignore 
the antiboson contribution to fluctuations. In addition, 
the effective charge shift is of order (P 2 / 'AI b )(8Mb / 'dfi), 
which is negligible by definition of nonrelativity, i.e., 
(P 2 /2M B <§C Mb)- This argument justifies the ordi- 
nary nonrclativistic treatment. Also at high tempcr- 

1 /3 

ature T <; -/V to ' t /Mb, the thermal fluctuation domi- 
nates over the quantum fluctuation because of the Bosc- 
enhancement in A^sr, i- e -; /s( a -') ~ T/u> for lo <C T, 
and we can safely ignore the quantum fluctuation. Thus, 
we expect that the effect of the quantum fluctuation be- 
comes less significant with going from the weak coupling 
BCS regime to the RBEC regime. 

In FIG. [151 we show how the (T c , G^ 1 )-relation is af- 
fected by the incorporation of the quantum fluctuation. 
We solved the number equation 



u3 



AW(^T) + N NS r{h,T) - 



(156) 

together with the Thouless condition Eq. ([55]) in (/i, T) to 
obtain the critical line (Hc(Gr),T c (Gb.)). The bold line 
is the result without the quantum fluctuation, while two 
thin lines show the results with two different choices for 
a cutoff Ab in fl q a; the lower (upper) line corresponds to 
Ab = 0.2A (0.3A). We conclude that the quantum fluctu- 
ation greatly affects the critical temperature in the BCS 
regime. This is a feature which is specific in the relativis- 
tic BCS state where the condition m ^> (|/i — m|, T c ) does 
not hold. As anticipated, the effect of quantum fluctua- 
tion becomes smaller for the (R)BEC regime where the 

1 /3 

critical temperature is relatively high (T c <; AT tot /Mb)- 
Let us finally explain why the quantum fluctuation in- 
creases the critical temperature in the BCS regime. The 
quark number from the vacuum fluctuation is negative 
as seen in Eq. (|155p . Then, to keep the quark number 
unchanged, the quark number from the mean field con- 
tribution Nyip([i,T) should increase. To realize this on 



27 




FIG. 15: The critical temperature versus the four-Fermi cou- 
pling. The bold curve shows the result without the quantum 
fluctuation, while two thin lines are the results when the quan- 
tum fluctuation is incorporated with two choices for cutoff, 
i.e., A B = 0.2A (down) and A s = 0.3A (up). 



the Thouless line in (/z, T)-plane, the line (fi, T c (/x)) along 
which the Thouless criterion is satisfied, the temperature 
should also increase because T c (/i) is an increasing func- 
tion as well as -/Vmf(a«,?c (/■*))■ So wc have to go higher \x 
(and T c (/jL)) in the (fx, T)-plane in order to keep the total 
quark number unchanged. 



VI. SUMMARY AND OUTLOOK 

We have made an extensive analysis on the static and 
dynamic aspects of the BCS/BEC crossover in a relativis- 
tic superfluid. We first developed the relativistic formula- 
tion of the Noziercs-Schimit-Rink framework for a Dirac 
fcrmion interacting with a four-Fermi point attraction. 
By carrying out the regularization of the dynamic pair 
susceptibility using the low energy expansion of the T- 
matrix in vacuum, we have performed a systematic exten- 
sion of the nonrelativistic Nozieres-Schimit-Rink (NSR) 
scheme [50| to a relativistic fermion system. Although 
we aimed particularly at a possible crossover in quark 
matter, our framework is general being not restricted 
to quark matter, and may be used for other relativistic 
fermion systems such as a possible neutrino superfluid- 
ity [74| inside compact stars at the early stage of their 
thermal evolution. We have found that in the relativistic 
case, there is an additional source of fluctuation which 
we have called the "quantum fluctuation" , and that our 
approach is consistent with the traditional Nozieres- 
Schmit-Rink framework when (m 3> |/J. — m\,T c ) where 
quantum fluctuation is absent. 

We have shown by the numerical calculation that three 



physically distinct regimes appear successively when the 
attraction is increased: In the weak coupling regime, the 
system is in the BCS phase where the critical temper- 
ature is exponentially smaller than the Fermi energy. 
When the coupling is increased beyond the unitary limit 
a s = ±00, the system gradually goes into the ordinary 
BEC phase where the increase of the critical temperature 
is little suppressed because of the appearance of bound 
states in the spectrum; in this regime, the increase in the 
attraction results mainly in the stabilization of boson. 
However, the critical temperature still shows a gentle in- 
crease due to the mass shift of in-medium bosons, which 
is highly in contrast to the nonrelativistic case. When the 
attractive coupling is increased further beyond the criti- 
cal coupling for the Majorana mass formation in vacuum, 
the system goes into the relativistic BEC (RBEC) phase 
[57[ | where the thermodynamics is dominated by anti- 
bosons, antifermions as well as fermions and bosons. To 
make physics of these three phases transparent, we dis- 
cussed how the spectral function, occupation numbers, 
Cooper pair size, and entropy density behave as a func- 
tion of attraction. 

We have also clarified the physical meaning of "quan- 
tum fluctuation" and have shown numerically that it 
greatly affects the critical temperature in the BCS 
regime, which is also in contrast to the nonrelativistic 
BCS phase. 

In addition to the study of the static aspects of the 
crossover, we have also studied the evolution of soft mode 
dynamics. By carrying out the low-energy and long- 
wavelength expansion of the dynamic pair susceptibility, 
we have seen how the effective theory for the dynamics of 
fluctuating pair field changes throughout the crossover. 
We saw that, in the weak coupling BCS regime, the dy- 
namics of soft modes is described by the time-dependent 
Ginzburg-Landau (TDGL) theory. Accordingly, the pair 
fluctuation is an ovcrdamped mode, which dissociates 
into two fermions in the Fermi sea [76|, [8(| ■ When the sys- 
tem goes into the BEC regime crossing the unitary limit, 
the TDGL theory is taken over by the Gross-Pitaevskii 
(GP) theory as in the nonrelativistic calculation (76[. In 
this regime, the low energy pair fluctuating field becomes 
a propagating mode which cannot decay into fermions 
because of the bound state gap. When the coupling is 
further increased and the system goes into the RBEC 
phase, the relativistic Gross-Pitaevskii (RGP) theory in 
turn describes a wide kincmatical region. In this region, 
the velocity of the fluctuation is large and the repulsive 
force between them is also large. 

Based on these low energy effective theories, we dis- 
cussed how the shear viscosity behaves throughout the 
BCS/BEC crossover. In the spirit that the short range 
quantum correlation is already taken into the low-energy 
coefficients in the effective theories, and also on the basis 
of the ansatz that soft modes are dilutely distributed in 
the phase space, we estimated the shear viscosity using 
the classical Kubo's formula and the Boltzmann equa- 
tion. Wc have found that the viscosity in the BCS regime 
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comes from the diffusion process of soft modes, while it 
comes from the binary scattering between soft modes in 
the (R)BEC regime. Our model calculation indicated 
that the soft mode contribution is important near the 
unitary limit; the shear viscosity to the entropy density 
ratio takes the minimum in the unitary regime, but the 
minimum bound of the viscosity proposed by Kovtun, 
Son and Starinets [TTJ] is kept owing to soft modes. 

Unfortunately, as the density is increased (ftp to), 
the Pauli-blocking effect against the formation of in- 
medium bosons, pushes the BCS/BEC crossover bound- 
ary to a larger attraction, and makes the crossover char- 
acteristics of the critical temperature somewhat ambigu- 
ous. Also a typical diquark coupling which is usually 
adopted in the literature of color superconductivities 
[33l.l34| is located well below the unitary point. Based on 
these points, we have concluded that the following "ex- 
otic" conditions are required to have the diquark BEC 
[55l . [HI in quark matter core in compact stars; (i) the 
attraction bet ween quarks is much larger than expected 
perturbatively [l02j , and (ii) the in- medium fermion mass 
is la rger than the perturbative hard dense loop calcula- 
tion [1031 ] and must be comparable to chemical potential. 

Investigating the transport properties for various can- 
didates of the color superconducting phase of qu ark m at- 



ter (104l | as well as the normal nuclear matter |105[, or 
exploring the BCS/BEC crossover with paying respect 
to more realistic conditions of quark matter remains an 
important issue. It may also be interesting to extend 
our framework so as to take into account the feedback of 
soft modes to the fermion propagation in medium, and 
to ask ho w th e fermion is modified near the critical tem- 
perature [l06| |. The self-consistent T-matrix theory [77| 
or the relativistic version of the Brueckner-Hartree-Fock 



theory [107| will be one of the best ways to study the 
in- medium fermion property affected by soft modes. 



Acknowledgments 

This work is motivated in part by my previous study 
[55l | performed in collaboration with Y. Nishida to whom 
I am very grateful. The first part of this paper is greatly 
indebted to fruitful discussions with him. I have also been 
greatly benefited from discussions with M. Kitazawa, 
T. Kunihiro and Y. Nemoto, which I acknowledge here. I 



am thankful to J. P. Blaizot, K. Itakura, E. Kolomcitscv, 
G. Nardulli, K. Ohnishi, M. Ruggieri, and K. Tsumura 
for helpful comments, discussions, and informing me of 
some references on this topic. I also thank T. Hatsuda 
for useful comments on pair size. The numerical calcula- 
tions were carried out on Altix3700 BX2 at YITP owned 
by Kyoto University. This work was supported in part 
by a Grant-in- Aid for the 21st Century COE "Center for 
Diversity and Universality in Physics" . 

APPENDIX A: IMAGINARY PART OF THE 
DYNAMIC PAIR SUSCEPTIBILITY 



The imaginary part of the dynamic pair susceptibility 
can be calculated as 



Im x m ,t(w-M<5, P) 



9[z> y/4m> + P*\ftf{z,P) 
-9[\z\<P]f^(z,P) 



(Al) 



[z<-y/4m? + P*\f^(z,P), 



where wc have defined a variable z = uj+2fi. We can show 
the identity fgf(z,P) = f-^' T (-z,P). The explicit 
forms for f^' T and f^ T can be found as 



v /z 2 -P 2 
4tt 



yjz 2 - P 2 - 4m 2 



2tvP 



-Tin 



l+e 



2-nP 



2ixP 



-Till 



l+e 



l + e 



(A2) 

/ld is present only for T ^ 0, and is called as the Landau 
damping, i.e., the Chcrenkov radiation (absorption) of 
soft mode from (to) the thermally excited quarks. We can 
easily check the identity Imx M ,T(0,P) = which states 
that the quasi-quark excitation on the Fermi surface is 
stable. 
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